C77-464/501 


and 

D.  L.  Mills 
Consultant 

University  of  California,  Irvine 
1 June  1977 


U.  S.  ARMY  RESEARCH  OFFICE 
ARO  Contract  No.  DAHC04-74-C-0024 
1 June  1974  to  1 June  1977 


■ApprrvverJ 

unlimited. 


for  public 


'■eJease ; 


distiibutfon 


Q- 

o 

CJ> 


an 

So 


Electronics  Research  Center 
Rockwell  International  Corporation 
3370  Miraloma  Avenue 
Anaheim,  California  92803 


•-J 


Unclassified 


security  classification  of  This  page  (Whtn  D»lm  Enitrtd) 


I 1/ 


I-  L - If  .9  - 


REPORT  DOCUMENTATION  PAGE 


\^C77-464/501 


2.  GOVT  ACCESSION  NO 


INVESTIGATION  OF^URFACE  Jlt’TICAL  WAVES  FOR  ; 
• OPTICAL  SIGNAL  PROCESSIN(T.  — t 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


T.yfJJif  report  a PERIOD  COVERED 

Final  technical  Report. 

1 Juno  1974  T June  1877 


7.  AUTHORf*; 

/ P 


J.  D./McMullen  0.  L.^Mills 


/< 


9 PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Electronics  Research  Center 

Rockwell  International  /6  ' 

3370  Miraloma  Ave.,  Anaheim,  CA  92803 


n.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Army  Research  Office  . 

Box  CM,  Duke  Station 

Durham,  North  Carolina  — 


/ / 


rA 


14.  monitoring  agency  name  » ADDRESSC// d///eren(  from  Ctmlrolling  Oftico) 


i PFRFOBminG  org.  report  NUMBET*- 

C/7-464/b01 


6.  CONTRACT  OR  GRANT  NUMBERfs; 


DAHC04-74-C-0024 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  & WORK  UNIT  NUMBERS 


D0161102B11B 


Ig  PFPAOT  nATLC 

^1  June  1977  / 


13.  NUMBER  OF  PAGES 


139 


15.  security  class,  (of  th/a  report) 

Unclassified 


15a.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  ST ATEMENT  r©/ f/j/*  ReporO 


r.f. 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (of  fhe  abafracf  er\tered  in  Block  20,  if  different  from  Report) 


10.  SUPPLEMENTARY  NOTES 


The  views  and  conclusions  contained  in  this  document  are  those  of  the 
authors,  and  should  not  be  interpreted  as  necessarily  representing  the  official 
policies,  either  expressed  or  implied,  of  the  Army  Research  Office  or  the 
United  States  Government. 


19.  KEY  WORDS  (Continue  on  reverse  aide  if  necessary  and  identify  by  block  number) 

Surface  Pol ari tons.  Surface  Optical  Waves,  Surface  Plasmons, 
CO2  Pulsed  Laser  Radar,  Frequency  Chirped  Pulse  Compression, 
Surface  Roughness,  Dispersive  Media,  Phase  Modulation. 


20.  abstract  (Continue  on  reverse  aide  If  necessary  and  identify  by  block  number) 


DD  , 3 1473  EDITION  OF  I NOV  65  IS  OBSOLETE 


Unclassified 


security  CLASSIFIC.'TION  of  this  page  flFhen 


1 


- / 


Data  Entered)  ^ 


security  CL*SSIFICATI0N  or  this  PACECWTufi  Dmim  Enfrmd) 


i 


ABSTRACT 

Propagation  and  temporal  compression  of  frequency-chi rped  COp  laser 
pulses  has  been  investigated,  wherein  a dispersive  optical  pulse  delay  line 
is  formed  usiig  dispersive  surface  and  bulk  phonon-polari ton  propagation 
modes  in  solids.  Criteria  are  developed  for  the  compression  of  optical 
pulses,  and  these  criteria  are  compared  to  the  relevant  group-dispersive 
properties  of  the  propagation  modes  to  determine  their  suitabilities  for 
performing  laser  pulse  compression.  Surface  phonon-polari ton  modes  on 
BeO:air  and  on  BeO:Ge  film:air  structures,  and  bulk  phonon-polaritons  in 
a -quartz  are  examined  as  dispersive  media  for  CO2  laser  pulses. 

Absorption  in  the  infrared-active  medium  has  been  shown  to  limit  the 
magnitude  of  group  dispersion  available,  in  addition  to  limiting  the 
attenuation  length  of  each  propagation  mode.  For  the  materials  considered, 
absorption  limits  the  application  of  this  approach  to  CO2  laser  pulses 
having  initial  widths  of  10  psec  or  smaller  with  initial  chirp  bandwidths 
in  the  THz  range,  if  attenuation  by  absorption  is  to  be  limited  to  50  dB. 

Both  the  narrow  pulse  width  and  large  chirp  bandwidth  requirements  preclude 
experimental  demonstration  with  present  CO2  laser  technology.. 

Consequently,  attention  has  been  focused  in  part  on  other  pulse 
compression  schemes  appropriate  for  CO2  lasers.  One  method  recognized 
as  promising  is  the  grating-pair  pulse  compression  filter  in  a highly 
dispersive  configuration.  An  exact  analytical  solution  has  been  developed 
for  chirped  pulses  compressed  by  this  type  of  strongly  dispersive  filter. 

An  alternative  scheme  has  been  explored,  using  active  modulation  of 
the  phase  velocity  of  a medium  to  compress  the  pulse  envelope. 

In  earlier  work  on  this  program  prism  and  grating  coupling  were 
investigated.  Attenuation  lengths  were  measured  for  surface  polaritons 
on  metals  and  BeO,  and  were  found  to  be  macroscopic.  Analysis  of  scattering  of 
surface  polaritons  by  surface  roughness  showed  that  acceptable  BeO  surface 
quality  could  be  produced  by  conventional  optics  polishing  technology. 
Interactions  between  optical  surface  waves  and  surface  elastic  waves  were 
investigated  as  a means  to  produce  the  desired  chirp  phase  modulation. 

The  influence  of  dielectric  overlayers  upon  propagation  characteristics 
of  surface  polaritons  was  studied,  as  means  to  modify  the  dispersion 
curves  for  particular  pulse  compression  applications. 

Nonlinear  optical  interactions  among  surface  polaritons  have  been  studied, 
including  experimental  attempts  to  observe  second  harmonic  generation  by  a 
10.6  pm  surface  plasmon  on  a metal  surface.  Recent  theoretical  work  (Ch.  VI) 
has  shown  that  the  second  harmonic  signal  to  be  expected  in  the  previously 
attempted  experiment  was  just  below  the  limit  of  detectability.  This 
experiment  is,  therefore,  being  attempted  again  with  appropriate  modifications, 
and  may  yield  the  first  observation  of  nonlinear  optical  interactions  among 
localized  surface  polaritons.  Positive  experimental  results,  if  forthcoming, 
will  be  reported. 
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I.  INTRODUCTION 

A.  Background 

The  purpose  of  this  program  is  to  investigate  those  properties  of 
surface  polaritons  which  may  be  applicable  for  signal  processing  at  CO^ 
laser  frequencies,  particularly  for  dispersive  optical  pulse  delay  lines  i 

to  compress  frequency  chirped  CO2  laser  pulses.  | 

In  this  application,  a pulse  having  just  the  appropriate  amount  of  I 

j 

of  frequency  chirp  becomes  compressed  after  propagating  through  a dispersive 
medium,  while  background  unchirped  optical  signals  become  dispersion- 
broadened  and  attenuated.  The  higher  peak  power  of  the  compressed  pulse, 
along  with  suppression  of  background  signals,  offers  improved  signal  ; 

detection,  while  the  narrow  width  of  the  compressed  pulse  gives  higher  | 

range  resolution.  ^ 

Conventional  methods  of  pulse  compression  for  pulsed  laser  radars 
require  converting  the  signal  to  electrical  and  acoustic  signals,  with 
large  iniertion  losses  or  bandwidth  limitations  and  sophisticated,  power-  ■ 

I 

consuming  associated  electronics.  Pulse  compression  performed  directly  ! 

at  optical  frequencies  eliminates  several  of  these  disadvantages.  The  ^ 

present  approach,  using  a dispersive  propagation  mode  in  a solid,  compares  j 

favorably  to  other  optical  techniques  such  as  grating  pairs  or  interferometers,  | 

because  of  the  relative  compactness  of  the  pulse  delay  line  and  the  possibility  j 

I 

of  using  ultra-wide  chirp  bandwidths  to  compress  broader  laser  pulses  than 
previously  possible.  Further,  surface  polariton  propagation  modes  in  over- 
layer film  structures  offer  the  possibility  for  constructing  a pair  of 
matched  filters,  so  that  both  chirp  and  compression  may  be  performed  by 
compact,  passive  elements  at  optical  frequencies. 


The  primary  problems  to  be  addressed  are  (1)  the  efficiency  with  which 
a chirped  laser  pulse  can  be  coupled  to  the  surface  propagation  mode,  and 
possible  distortion  of  the  phase  modulation  (chirp)  and  pulse  shape  during 
the  coupling  process,  (2)  propagation  losses  arising  from  absorption  and 
from  roughness-induced  scattering  of  the  surface  polariton,  and  (3)  the 
amount  of  group  velocity  dispersion  attainable  for  surface  propagation 
modes  at  CO^  laser  frequencies. 

Significant  Results 

Both  prism  and  grating  couplers  were  investigated  as  means  to  couple  a 
CO2  laser  beam  to  surface  propagation  modes  on  BeO,  quartz  and  metal  sub- 
strates. Coupling  efficiencies  of  ten  percent  were  measured  for  prism 
couplers  used  to  launch  surface  plasmons  on  metals,  and  about  one  percent 
for  generation  of  a surface  phonon-polariton  at  10.6  ym  on  BeO.  Attenuation 
lengths  at  10.6  ym  were  typically  one  to  8 mm  for  surface  plasmons  on  stain- 
less steel,  copper  and  gold,  and  120  ym  for  the  ordinary  surface  phonon- 
polariton  on  BeO. 

Methods  to  modulate  the  phase  of  a surface  optical  wave  have  been 
investigated,  including  both  electro-optic  modulation  and  interactions  with 
Rayleigh  waves.  These  interactions  offer  the  possibility  to  chirp  the 
pulse  while  it  is  propagating  as  a surface  wave. 

Nonlinear  optical  interactions  among  surface  polaritons  have  been 
studied.  Experiments  to  generate  a second  harmonic  surface  plasmon  by  a 
fundamental  surface  plasmon  at  10.6  ym  did  not  yield  a detectable  harmonic 
signal . 

The  influence  of  dielectric  overlayers  upon  the  dispersion  curves  of 
surface  phonon-polaritons  has  been  examined  in  detail  as  a means  to  modify 
the  curvature  of  the  dispersion  curves  of  surface  phonon-polaritons. 
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This  approach  permits  the  design  of  a matched-filter  pair  for  expansion  and 
compression  of  laser  pulses. 


The  theory  for  chirped  pulse  compression  has  been  generalized  to 
include  strongly  frequency-dependent  absorption  in  the  dispersive  medium. 

The  role  of  absorption  has  been  quantitatively  evaluated  for  both  surface 
and  bulk  phonon-polaritons  at  CO2  laser  frequencies.  It  has  been  shown 
that  absorption  limits  both  the  attenuation  length  of  the  mode  and  the 
magnitude  of  group  velocity  dispersion  attainable.  Because  of  absorption, 
pulse  compression  using  surface  phonon-polaritons  on  BeO  or  bulk  phonon- 
polaritons  in  quartz,  requires  that  the  initial  CO^  laser  pulse  be  chirped 
over  a THz  range,  or  an  entire  P or  R branch,  and  that  even  then  the  initial 
jiulse  width  be  20  psec  or  less.  These  parameters  preclude  experimental 
demonstration  with  existing  COp  laser  technology.  These  results  are  expanded 
upon  in  Section  II  of  this  report. 

^ The  theory  for  compression  of  a chirped  pulse  in  a strongly  dispersive 
mediurti  has  been  examined,  and  new,  closed-form  expressions  developed  for  the 
shape  of  the  compressed  pulse.  In  this  treatment  it  is  assumed  that  both 
the  group  velocity  v^^  and  its  first  frequency  derivative  dv^^/doj vary 
significantly  with  frequency  over  the  bandwidth  of  the  incident  chirped 
pulse.  An  approach  has  been  identified  for  compression  of  CO2  laser  pulses 
using  a grating-pair  filter  in  a highly  dispersive  configuration.  Elements 
of  this  analysis  follow  in  Section  III  of  this  report. 

Active  modulation  of  the  refractive  index  of  a crystal  has  been  examined 
as  a means  to  compress  optical  pulses,  as  outlined  in  Section  V . Criteria 

are  developed  for  active  phase  modulation  and  for  distortion  of  the  pulse 
envelope  upon  passing  through  a modulator  of  finite  thickness. 
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C.  List  of  Publications  and  Technical  Reports 

The  following  papers  have  been  published  or  submitted  for  publication 

as  results  of  this  program: 

D.  L.  Mills,  "Theory  of  Second  Harmonic  Generation  by  Surface  Pol ari tons 
on  Metals, Sol  id  State  Communic.  (to  be  published). 

J.  D.  McMullen,  "Chirped  Pulse  Compression  in  Strongly  Dispersive 
Media,"  J.  Appl.  Phys.  (to  be  published). 

J.  D.  McMullen,  "Analysis  of  Frequency-Chirped  Optical  Pulses  using 
Dispersive  Bulk  and  Surface  Phonon-Polariton  Propagation  Modes 
at  CO^  Laser  Frequencies,"  J.  Appl.  Phys.  (manuscript  No.  1209R, 
to  be  published). 

D.  L.  Mills,  "Theory  of  Surface  Polariton  Generation  through  Grating 

Couplers,"  J.  Appl.  Phys.  (manuscript  No.  1210R,  to  be  published). 

D.  L.  Mills,  "The  Interaction  of  Surface  Polaritons  with  Periodic 

Surface  Structures;  Rayleigh  Waves  and  Gratings,"  Phys.  Rev.  B 
1^,  3097  (1977). 

J.  D.  McMullen,  "Recursive  Solution  of  the  Dispersion  Relation  of  a 

Surface  Polariton  in  a Layered  Structure,"  J.  Opt.  Soc.  Am  698 
(1977). 

J.  D.  McMullen,  "Perturbation  Analysis  of  the  Influence  of  a Thin 

Dielectric  Film  on  the  Propagation  Constants  of  a Surface  Optical 
Wave,"  J.  Opt.  Soc.  Am.  6^,  1253  (1975). 

D.  L.  Mills,  "The  Attenuation  of  Surface  Polaritons  by  Surface  Roughness," 
Phys.  Rev.  B 12,  4036  (1975). 

J.  D.  McMullen,  "Propagation  Length  of  Surface  Polaritons  at  10.6  pm 

on  BeO  and  Metal  Surfaces,"  Solid  State  Communic.  17,  331  (1975). 

J.  D.  McMullen,  "Surface  Polariton  Attenuation  Length  Measurements  at 
10.6  ym  by  Direct  Excitation  and  Decoupling  and  by  the  Total 
Internal  Reflection  Method,"  Bull.  Am.  Phys.  Soc.  20^,  45  (1975). 

The  following  Technical  Reports,  entitled  "Surface  Electromagnetic  Waves 

for  Optical  Signal  Processing,"  have  been  submitted  to  A.R.O.  during  this 

program: 

Final  Technical  Report,  Period  1 June  1974  to  1 June  1977,  Report  No. 
C77-464/501  (1  June  1977). 
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Fifth  Semiannual  Progress  Report,  Period  30  June  1976  to  31  December 
1976,  Report  No.  C76-700. 5-501  (31  Dec.  1976). 

Fourth  Semiannual  Progress  Report,  Period  31  December  1975  to  30  June 
1976,  Report  No.  C76-700-501  (30  June  1976). 

Follow-on  Proposal,  No.  T76-275/501  (16  March  1976). 

Third  Semiannual  Progress  Report,  Period  1 June  1974  to  31  May  1976, 

Report  No.  C74-930. 5/501  (30  September  1975). 

Final  Technical  Report,  Period  1 June  1974  to  31  May  1975,  Report  No. 
C74-930. 4/501  (31  May  1975). 

Third  Letter  Progress  Report,  Period  15  December  1974  to  15  March  1975, 
Report  No.  C74-930. 3/501  (15  March  1975). 

Second  Quarterly  Letter  Technical  Report,  Period  15  September  1974 

to  15  December  1974,  Report  No.  C74-930. 2/501  (15  December  1974). 
Progress  Report  No.  1,  Period  1 June  1974  to  30  September  1974,  Report 
No.  C73-930. 1/501  (15  September  1974). 

D.  Abstracts  of  Publications 

Abstracts  of  papers  published  during  this  contract  period  are 
listed  as  follows.  Complete  copies  of  papers  submitted,  but  not  as 
yet  published,  follow  as  sections  II  through  VI of  this  report.  These 
sections  contain  new  results  obtained  during  the  present  reporting  period. 


■ 


I 
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Interaction  of  surface  polaritons  with  periodic  surface  structures; 

Rayleigh  waves  and  gratings^ 

D.  L.  Mills 
Department  of  Physics 
University  of  California 
Irvine,  California  92664 

ABSTRACT 

We  describe  the  interaction  of  surface  polaritons  with  Rayleigh 
surface  waves  for  the  case  of  arbitrary  angle  between  the  surface- 
polariton  propagation  direction  and  that  of  the  Rayleigh  wave.  We 
presume  the  surface  polariton  couples  to  the  Rayleigh  wave  through  the 
periodic  displacement  of  the  surface  produced  by  it,  and  also  through 
the  modulation  of  the  dielectric  constant  of  the  substrate  by  the 
elasto-optic  effect.  We  take  the  unstressed  substrate  to  be  an  isotropic 
dielectric,  and  no  special  assumption  about  its  elastic  constants  is 
required.  The  results  here  may  be  applied  also  to  the  interaction  of 
surface  polaritons  witi'  periodic  gratings,  provided  the  depth  of  the 
grooves  is  small.  We  obtain  expressions  for  the  amplitude  of  the 
diffracted  electromagnetic  waves  produced  by  interaction  of  the  surface 
polariton  with  a grating  or  Rayleigh  wave.  The  diffracted  beams  may 
describe  either  fields  localized  near  the  surface  or  radiation  induced 
by  interaction  with  the  periodic  distrubance,  depending  on  the  kinematics 
of  the  interaction  process.  A small  modification  of  these  formulas  may 
be  used  to  compute  the  amplitude  of  a surface-polari ton  wave  generated 
by  interaction  of  incoming  electromagnetic  radiation  with  a periodic 
structure.  Finally,  we  obtain  a remarkably  simple  analytic  expression 
for  the  surface-polariton  dispersion  relation  in  the  near  vicinity 
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of  gaps  opened  up  by  interaction  with  the  periodic  surface  structure. 

This  expression  includes  the  effect  of  elasto-optic  coupling  and  may 
be  applied  to  propagation  at  an  arbitrary  angle  to  a set  of  grating 
lines,  or  to  the  propagation  direction  of  a Rayleigh  wave.  The  exper- 
imental study  of  these  gaps  offers  the  possibility  of  measuring  the 
magnitude  and  frequency  variation  of  the  elements  of  the  photoelastic 
tensor  of  the  substrate  in  situations  where  these  quantities  are  difficult 
to  measure  by  conventional  methods;  i.e.,  for  metals  at  infrared  frequencies. 


*Research  supported  by  Contract  No.  DAHC04-74-C-0024,  U.  S.  Army  Research 
Office,  Durham,  N.  C. 
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Recursive  solution  of  the  dispersion  relation 
of  a surface  polariton  in  a layered  structure.* 


J.  D.  McMullen 

Electronics  Research  Division 
Rockwell  International  Corporation 
Anaheim,  California  92803 


Abstract 

The  transcendental  equation  for  the  propagation  constants  of  a surface 
polariton  in  a layered  structure  is  written  in  a form  which  is  particularly 
amenable  to  solution  by  a recursive  method.  A numerical  example  is  given 
for  the  surface  phonon-polariton  on  a structure  consisting  of  a Ge  film  on 
BeO. 

Index  Headings:  Surface  optical  waves;  films. 


*This  research  was  supported  by  the  U.S.  Army  Research  Office, 
under  Contract  DAHC04-74-C-0024. 
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PERTURBATION  ANALYSIS  OF  THE  INFLUENCE  OF  A THIN  DIELECTRIC  FILM 
ON  THE  PROPAGATION  CONSTANTS  OF  A SURFACE  OPTICAL  WAVE* 

J.  D.  McMullen 

Rockwell  International 
Anaheim,  California  92803 

ABSTRACT 

Simplified  analytic  expressions  are  obtained  for  the  phase  velocity  and 
propagation  length  of  a surface  polariton  on  a layered  structure  that 
consists  of  a dielectric  film  on  a surface-active  material. 

Index  Headings:  Surface-guided  waves;  Films;  Integrated  optics. 


*This  research  was  supported  by  the  U.  S.  Army  Research  Office,  under 
Contract  No.  DAHC04-74-C-0024 


THE  ATTENUATION  OF  SURFACE  POLARITONS 


BY  SURFACE  ROUGHNESS‘S 


D.  L.  Mills 
Department  of  Physics 
University  of  California 
Irvine,  California  92664 

ABSTRACT 

This  paper  presents  a theoretical  description  of  the  attenuation 
of  surface  polaritons  by  roughness  on  the  surface.  In  the  presence 
of  surface  roughness,  and  in  a frequency  region  where  the  dielectric 
constant  is  negative,  the  surface  polariton  is  attenuated  by  two 
processes.  It  may  lose  energy  by  radiating  into  the  vacuum  or  by 
scattering  into  other  surface  polariton  states.  Through  application 
of  a formalism  developed  recently  to  describe  roughness  induced 
scattering  and  absorption  of  a plane  electromagnetic  wave  incident 
on  a surface,  we  obtain  expressions  for  the  contribution  to  the 
attenuation  rate  of  the  surface  polariton  from  the  two  processes 
described  above.  We  examine  the  relative  importance  of  the  two 
processes  for  surface  polaritons  on  semiconductor  surfaces,  and 
on  a nearly  free  electron  metal  at  infrared  frequencies. 


tResearch  supported  by  Contract  DAHC04-74-C-0024,  U.S,  Army 
Research  Office,  Durham,  North  Carolina 
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PROPAGATION  LENGTHS  OF  SURFACE  POLARITONS  AT 
10.6  MICRONS  ON  BeO  AND  METAL  SUBSTRATES* 

J.  D.  McMullen 

Electronics  Research  Division 
Rockwell  International  Corporation 
Anaheim,  Ca  92803,  U.  S.  A. 

ABSTRACT 

Propagation  lengths  at  10.6  ym  have  been  measured  to  be  120  ym  for  the 
surface  phonon-polariton  on  BeO  and  5.0,  7.4  and  1.6  mm  for  surface 
plasmons  on  gold,  copper  and  stainless  steel.  A CO2  laser  beam  was  coupled 
to  the  surface  mode,  propagated  along  the  surface,  and  decoupled  using 
the  two-prism  method. 
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SURFACE  POLARITON  ATTENUATION  LENGTH  MEASUREMENTS  AT  10.6  MICRONS  BY  DIRECT 
EXCITATION  AND  DECOUPLING  AND  BY  THE  TOTAL  INTERNAL  REFLECTION  METHOD 

J.  D.  McMullen 

Rockwell  International,  Anaheim,  CA 
ABSTRACT 

Propagation  lengths  were  measured  for  surface  polaritons  by  two  methods. 

CO^  laser  radiation  at  10.6  microns  was  used  to  directly  excite  surface 

polaritons  by  total  internal  reflection  from  an  adjacent  dielectric  prism. ^ 

For  one  type  of  experiment,  damping  of  the  surface  polariton  is  deduced  from 

the  angular  width  of  the  coupling  angle  for  which  the  totally  reflected 

intensity  is  minimum.  In  the  other  measurements,  the  surface  polariton  is 

launched  by  coupling  from  a right-angle  prism  and  decoupled  by  another  orism 

2 

in  the  manner  introduced  by  Schoenwald,  et  al.  Results  of  the  two  methods 
are  compared  for  surface  plasmons  on  metals  and  the  ordinary  surface  phonon- 
polariton  on  BeO. 
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II.  ANALYSIS  OF  COMPRESSION  OF  FREQUENCY-CHIRPED  OPTICAL  PULSES 
USING  DISPERSIVE  BULk  AND  SURFACE  PHONON-POLARITON 
PROPAGATION  MODES  AT  CO2  LASER  FREQUENCIES* 


J.  D.  McMullen 
Electronics  Research  Center 
Rockwell  International  Corporation 
Anaheim,  California  92803 


Abstract 

An  analysis  is  given  for  the  compression  of  linearly  chirped  optical 
pulses  propagating  in  dispersive,  absorbing  media.  Numerical  calculations 
are  reported  for  frequency-chirped  CO2  laser  pulses  which  propagate  as 
phonon-polaritons  in  a -quartz  or  as  surface  phonon-polari tons  on  BeO, 
including  the  influence  of  thin  dielectric  films  on  the  surface  propagation 
mode.  Quantitative  evaluation  is  made  of  the  use  of  these  modes  to  form 
dispersive  optical  pulse  delay  lines.  It  is  shown  that  absorption  limits 
both  the  attenuation  length  and  the  curvature  of  the  dispersion  curves. 

For  the  materials  considered,  this  approach  is  restricted  to  10  picosecond 
pulses  having  THz  chirp  bandwidths,  with  large  accompanying  absorption 
losses. 


*This  research  was  supported  by  the  U.  S.  Army  Research  Office 
under  Contract  No.  DAHC04-74-C-0024 


It  is  well  recognized  that  a pulse  may  be  temporally  compressed  by 

appropriately  modulating  its  carrier  frequency,  then  propagating  it  through 
1 -4 

a dispersive  medium.  The  purpose  of  this  paper  is  to  evaluate  the  suit- 
ability of  using  surface  and  bulk  phonon-polari ton  propagation  modes  in 
crystals  to  form  dispersive  pulse  delay  lines  at  infrared  frequencies. 
Particular  attention  is  given  to  the  possibility  of  compressing  CO2  laser 
pulses  at  frequencies  between  900  cm"^  and  1100  cm'\ 

Phonon-polaritons  are  the  electromagnetic  modes  of  propagation  in  a 

crystal  at  frequencies  close  to  the  transverse  optic  (TO)  phonon  frequency 
5 

of  the  lattice.  The  dispersion  curve  for  the  coupled  excitation,  or 
phonon-polari ton,  bends  asymptotically  toward  large  values  of  wave  vector 
for  frequencies  increasing  toward  the  propagation  stopband,  because  of 
near- resonant  coupling  between  the  optical  electric  field  and  the  infrared 
active  phonon.  A surface  optical  wave  may  propagate  nonradiati vely  along 
the  surface  of  the  infrared  active  crystal  at  frequencies  just  above 
the  TO  phonon  frequency.^ The  dispersion  curve  for  the  surface  phonon- 
polari  ton  also  bends  toward  large  values  of  wave  vector,  with  a limiting 
frequency  just  below  the  longitudinal  optic  (LO)  phonon  frequency  of  the 
crystal.  The  bending  of  the  dispersion  curves  indicates  that  group  velo- 
city varies  strongly  with  frequency  over  these  ranges  of  frequency. 
Therefore,  different  frequency  components  within  the  spectrum  of  an  incident 
pulse  are  delayed  by  different  amounts  as  the  pulse  propagates  through 
the  dispersive  medium.  The  output  pulse  envelope  is  then  distorted  and, 
under  certain  conditions,  may  be  temporally  compressed. 
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COMPRESSION  OF  LINEARLY  CHIRPED  PULSES  IN  DISPERSIVE,  ABSORBING  MEDIA 

Compression  of  frequency  chirped  pulses  is  used  to  improve  the 

detection  sensitivity  and  range  resolution  of  pulsed  radarJ'^  Recently, 

emphasis  has  been  placed  upon  extending  pulse  compression  methods  to 

optical  frequencies  for  laser  communications  and  radar.  ' Frequency 

modulated  laser  pulses  suitable  for  compression  have  been  produced  by 

electro-optic  phase  modulation,  ’ and  by  self-phase  modulation  of  a 

pulse  propagating  in  a nonlinear  medium. Compression  of  frequency 

chirped  laser  pulses  has  been  investigated  using  dispersion  of  the 

refractive  index  of  a liquid  ^ or  of  an  alkali  vapor. Pulse 

compression  filters  have  also  been  demonstrated  using  diffraction  grat- 
is 19-<I2  23  24 

mgs  ’ or  a Gires-Tournoi s interferometer.  ’ Other  methods 

for  obtaining  compressed  laser  pulses  have  included  active  modulation  of 

pc  pc 

the  refractive  index  and  its  dispersion  using  a Stark  cell,  ’ and  an 

27 

intra-cavity  saturable  absorber  in  a mode-locked  CO2  laser. 

In  this  paper  the  propagation  of  a frequency-chi rped  optical  pulse 
through  an  absorbing,  dispersive  medium  is  discussed  within  the  frame- 
work of  linear  dispersion  theory.  The  theoretical  approach  is 

28-  30 

similar  to  that  given  previously  for  a non-chirped  incident  pulse, 

1-3  9 12 

and  for  dispersive  but  non-absorbing  media.  ’ ’ The  present  treat- 
ment is  generalized  to  include  the  possibilities  of  a non-optimi zed 
initial  chirp,  a finite  attenuation  length  in  the  absorbing  medium,  and 

the  temporal  shaping  influence  of  strongly  frequency-dependent  absorp- 

31  -33 

tion  as  well  as  dispersion.  Criteria  are  established  for  the 
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Initial  chirp  rate  required  to  obtain  a given  pulse-width  compression 
ratio.  In  terms  of  the  incident  pulse  width  and  the  frequency  dependence 
of  the  complex  wave  vector  k(to).  Quantitative  evaluation  is  then  given 
for  chirped  CO2  laser  pulses  propagating  as  bulk  phonon-polaritons  in 
a-quartz  or  as  surface  phonon-polaritons  on  BeO. 

Consider  an  optical  pulse  having  electric  field  E(x,t)  which  propa- 
gates in  the  +x  direction  in  a linear,  dispersive  medium.  The  linearly 
chirped  incident  pulse,  specified  at  a point  x=0  just  inside  the  boundary 
of  the  dispersive  medium,  is  of  the  form: 


E(0,t)  = Eq  exp 


cos 


M t + 6(jJ_ 

0 m 


0) 


where  (rad. /sec)  is  the  average  optical  carrier  frequency,  and  (rad. /sec) 
is  the  applied  range  of  linear  frequency  modulation  (chirp)  over  the  full 
width  T of  the  intensity  envelope.  The  assumption  of  a Gaussian  envelope 
greatly  simplifies  the  analysis,  permitting  a closed-form  solution  for  the 
output  field  E(x,t)  while  retaining  the  important  aspects  of  chirped  pulse 
propagation  and  compression.  The  frequency  spectrum  ^(0,u)  of  the  incident 
pulse,  given  by  the  Fourier  transform  of  E(0,t),  is: 


^(0,(jj) 


exp 


-2(aj-a)^)^ 


exp 


-2(to+o)  ) 

0 

5 o).(5  0). 

0 0 


— j — V + 1 i tan 
- 1 ) 2 
m' 


(2) 
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( 


where 


1/2 


(3) 


is  the  full  width  at  e"^  of  peak  for  the  power  spectrum  |^(0,oj)|^,  and 
<5^0=  4/T  is  the  full  width  at  e"^  of  peak  for  the  power  spectrum  of 
the  unmodulated  pulse  envelope. 

According  to  linear  dispersion  theory,  the  electric  field  E(x,t)  after 
propagating  a distance  x through  the  dispersive  medium  is  given  by; 


/w 

dci)  1^(0, cj)  exp[+ik(cj)x  - iwt] 


(4) 


It  is  assumed  that  ^(0,w)  is  localized  to  frequencies  in  the  vicinities  of 
oj^+uJq  and  In  order  to  ensure  that  the  pulse  amplitude  diminishes 

to  zero  in  the  limit  x-h-®  in  an  absorbing  medium,  the  positive  sign  for 
+k(oj)  must  be  chosen  when  performing  the  integral  of  Eq.  (4)  over  the 
term  in  ^!(0,uj)  peaked  about  'jo^+uj^  and  -k(cu)  is  chosen  for  the  term  in 
1^(0, oj)  peaked  about  For  most  cases  of  interest  the  variation  of 

k(cj)  with  frequency  is  sufficiently  weak  that  k(<jo)  may  be  expanded  in  a 
Taylor  series  of  the  form 


k(u)x  = kgX  + 


(5a) 


about  u=+(jo  , and 
0 


k(w)x  = k^x  - T^(co+o)^)  + I T2^(ja+u)Q)^  + 


(5b) 


! 


J 
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( 


about  where  k^=k(uj^)  and  the  complex  parameters  i-j  and  i2  are 

defined  by; 

X , (5c) 

0 

and 

X . (5d) 

U)=+UJ 

0 

In  the  absence  of  absorption,  t-|  is  simply  the  group  delay  time  for  the 
pulse  envelope,  while  T2,  defined  here  as  the  "quadratic  dispersion  para- 
meter," characterizes  the  relative  delay  time  among  frequency  components 
after  propagation  through  the  dispersive  medium.  At  frequencies  for  which 
Im[k(u)]  varies  rapidly  with  frequency,  t-j  and  T2  are  more  generally  com- 
plex constants  which  are  related  to  the  slope  and  curvature  of  the  curves 
for  Re[k(o))]  and  Im[k((i))],  respectively. 

For  many  cases  of  interest  the  series  expansions  for  k(ijj)  may  be 
truncated  after  the  second-order  term.  If  the  frequency  dependence  of 
Im[k(u)]  is  negligibly  small,  then  the  electric  field  waveform  at  position 
X is  evaluated  to  be: 


(x^)^  = 4 d^k 


T - 

^1  = -j: 


^ ^ 
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In  the  absence  of  strongly  frequency-dependent  absorption,  so  that 
Im[k(uj)]  = Im[k^],  the  envelope  and  phase  portions  of  E(x,t)  in  Eq.  (6) 
are  readily  distinguishable.  E(x,t)  has  a Gaussian  envelope  of  full 

width  at  e’^  of  peak  intensity  given  by:  ! 


X 


I + m 


2 


+ ("2 


(7) 


in  agreement  with  Eq.  (22)  of  Reference  12.  The  peak  intensity  of  the 
pulse  at  position  x decreases  or  increases  by  the  factor  7/T^  as  the  pulse 
becomes  either  dispersion-broadened  or  temporally  compressed,  respectively. 
Additionally,  the  pulse  intensity  is  reduced  by  a factor  exp( -2Im(k^)x) 
because  of  absorption  by  the  medium. 

It  is  apparent  from  Eq.  (7)  that  if  no  frequency  chirp  is  applied  to 
the  incident  pulse  (6aj^=0),  the  output  pulse  E(x,t)  becomes  dispersion- 
broadened,  as  illustrated  by  the  dashed  curves  in  Figure  1.  The  phase  of 
E(x,t)  then  contains  a linear  frequency  chirp  covering  a range  of  frequen- 
cies ‘5*^0  width  T^  of  the  output  pulse  intensity 

2 

envelope.  However,  when  a frequency  chirp  <5w^=-6uj^(T/t2)  is  applied  to 

the  incident  pulse,  the  output  pulse  becomes  compressed  to  a minimum  pos- 

1 3 12 

sible  width  upon  propagation  through  the  dispersive  medium,  ’ ’ as 
illustrated  by  the  dotted  limiting  curve  in  Fig.  1.  The  width  of  the 
optimally  compressed  pulse  is  T^=4/l'5w^|,  and  the  phase  of  E(x,t)  contains 
a linear  frequency  chirp  during  full  width  T^  over  a range  of  frequencies 
5(jJq.  The  rate  of  chirp  of  the  optimally  compressed  pulse  is  equal  in  mag- 
nitude but  opposite  in  sign  to  the  rate  of  chirp  for  the  incident  pulse. 
The  intermediate  cases  forju^<0  represented  in  Fig.  1 by  solid  curves 
show  that  large  compression  ratios  can  be  obtained  for  sufficiently  large 
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chirp  bandwidths  and  sufficiently  small  quadratic  dispersion  parameter 
However,  optimum  compression  requires  an  increasingly  delicate  balance 
among  the  parameters  60)^^^  and  T of  the  incident  pulse  and  disper- 

sive medium  if  large  compression  ratios  are  to  be  attained. 

When  the  average  carrier  frequency  is  very  close  to  the  propaga- 
tion stopband  for  either  a bulk  or  surface  phonon-polariton,  the  frequency 
dependence  of  Im[k(u))]  may  also  contribute  significantly  to  pulse  shaping. 
The  output  pulse  width  is  then  given  more  generally  by: 


Re{T2^)l 


5u) 


^ ^ " [siZ  ~7T- 


(8) 


The  minimum  possible  output  pulse  width  is  then  given  by: 


CsJ 

X 

h- 

60J 

m 

^“o'' 

r 

6(jj 

L J 

T 

L 0 

m j 

when  the  input  chirp  is  adjusted  to  an  optimum  value  of: 

^2 


(5oj_ 


1 

opt. 


Re(t2^) 


1 + 


Im(T2^) 


(9b) 


This  reduces  to  the  previously  discussed  result  for  optimum  compression  in 
the  limit  that  frequency-dependent  absorption  becomes  sufficiently  weak, 
or  Im[k{.j)]  = Im(k^) . 
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When  Im(t2  ) becomes  comparable  to  T , as  may  occur  when  is  just  below 
the  propagation  stopband,  then  the  chirp  bandwidth  required  for  optimum 
compression  is  larger  and  the  attainable  pulse  compression  ratio  is 
reduced  because  of  frequency-dependent  absorption. 

BULK  PHONON-POLARITON  MODES  IN  QUARTZ 

The  preceding  discussion  has  shown  how  the  curvatures  of  the  curves 

for  Re[k((jj)]  and  Im[k(o))]  are  related  to  temporal  shaping  of  the  pulse 

envelope.  For  example.  Fig.  1 illustrates  that  the  quadratic  dispersion 

parameter  must  be  comparable  to  the  initial  pulse  width  if  a non-chirped 

pulse  is  to  be  significantly  dispersion-broadened.  Also,  when  the  chirp 

rate  in  the  incident  pulse  is  optimized  Eq.  (7)  shows  the  minimum  pulse 

2 

width  obtainable  is  /T.  Therefore,  the  quadratic  dispersion  parameter 
i2  must  be  calculated  for  bulk  and  surface  phonon-pol ari ton  propagation 
modes  in  order  to  evaluate  their  potential  use  for  optical  pulse  compres- 
sion. 

Consider  first  the  ordinary  phonon-pol ari ton  modes  in  a-quartz.  The 
dispersion  curves  illustrated  in  Fig.  2 are  given  by: 

k(-^)  = ^ (TO) 

where  the  complex,  frequency-dependent  dielectric  constant  e(oj)  for  the 
ordinary  ray  is  of  the  form: 


N c , 2 


= eH  + 53  — rV 

m=l  ^‘"m 


(11) 
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Sj^,  and  are  the  oscillator  strength,  transverse  optic  phonon  frequency, 

and  anharmonic  damping  rate,  respectively,  for  the  m'th  of  N ordinary 

phonon-polariton  branches,  and  e(«>)  is  the  ordinary  optical  dielectric 

constant  at  frequencies  well  above  the  lattice  resonance  but  still  well 

below  electronic  bandgap  transitions.  Data  for  the  parameters  used  to 

calculate  e(u)  are  taken  from  infrared  reflectivity  measurements  by 

34 

Spitzer  and  Kleinman.  The  solid  curves  of  Fig.  Z represent  phonon- 
polariton  dispersion  in  the  limit  of  no  damping  (y,^=0j  ni=l . 2, ...7), 
while  the  dashed  curve  for  Re[k(w)]  and  the  dotted  curve  for  Im[k(ijj)] 
include  the  influence  of  anharmonic  damping  at  300°K.  The  undamped  case 
is  presented  as  a limiting  case  of  theoretical  interest,  since  the  strong 
bending  of  the  undamped  dispersion  curve  stimulated  the  present  study. 
However,  an  assessment  of  practical  applications  clearly  must  include  the 
influence  of  absorption. 

The  parameter  T2  was  calculated  by  computer,  using  analytic  expres- 

2 2 

sions  for  the  derivatives  dk/doj  and  d k/doj  derived  from  Eqs.  10  and  11. 

These  expressions  are  somewhat  lengthy,  and  are  therefore  not  reproduced 

here.  Figure  3 illustrates  the  variation  of  with  average  carrier 

frequency  for  an  assumed  propagation  distance  of  x=lcm.  For  a smaller 

propagation  distance  x,  as  may  be  required  by  a small  attenuation 

length,  the  value  of  the  ordinate  must  be  multiplied  by  a factor 
1 /2 

[x(cm)]  ' . The  solid  curve  of  Fig.  3 represents  the  limit  of  no 

2 1/2 

damping.  With  absorption  included,  the  quantities  |Re(r2  )|  and 
2 1 /2 

|Im(T2  )|  ' are  illustrated  by  dashed  and  dotted  curves,  respec- 

2 1/2 

tively.  The  dashed  curve  for  |Re(T2  )|  follows  closely  the  solid 


r 
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I 

I 
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1075  cm  in  addition  to  a peak  for  negative  values  of  Re(t2  ) at 

1071  cm  . Negative  values  for  Re(i2  ) occur  between  1067  and 

1073  cm  ^ and  for  a range  of  frequencies  below  889  cm'\  for  which  the 

Re[k(uj)]  curve  in  Fig.  2 is  concave  toward  decreasing  values  of  k. 

p 1 

Maxima  occur  for  positive  values  of  Im(x2  ) at  1069  and  1080  cm"  and  a peak 

for  negative  values  at  1073  cm”^  frequency. 

As  an  example,  consider  the  propagation  of  a chirped  laser  pulse 

having  an  average  carrier  frequency  of  925  cm'\  corresponding  to  the 

(00°1 ) -( 10°0)  P(40)  emission  line  of  a CO2  laser.  The  attenuation  length 

[2Im(k)]  at  this  frequency  is  46  qm,  while,  from  Fig.  3,  j Re(x2  ) | '^  = 

5.0x10  sec  and  | Im(x2  )|  ' = 2.7x10  sec  for  x=lcm.  A quartz  crystal 

0.5  mm  thick  then  gives  j Re(x2^)  | = 1.1x10'^^  sec  and  | Im(x2^)  I = 

-12 

0.6x10  sec,  along  with  a reduction  in  peak  intensity  by  a factor 

exp[-21m(k^)x]  = 2x10  due  to  absorption.  It  is  apparent  from  the 

magnitude  of  X2  that  a pulse  which  is  to  be  optimally  compressed  in  this 

medium  should  have  an  initial  width  of  a few  picoseconds  and  requires  a 

chirp  bandwidth  of  the  order  of  a few  THz.  A chirp  bandwidth  of  this 

magnitude  would  require  a high-pressure  CO2  laser,  for  which  the  emission 

lines  of  the  (00°1 ) -( 10°0)  band  P branch  are  broadened  sufficiently  to 

35-39 

merge  into  a continuous  gain  spectrum.  For  a chirp  range  of  -1.5  THz, 

a 2.9  psec  pulse  can  be  compressed  to  0.43  psec  in  width,  with  a gain  factor 
T/T^  = 6.7  in  peak  intensity.  If  a broader  incident  pulse  is  used,  so 
that  conditions  for  pulse  compression  are  no  longer  optimized,  the  same 
chirp  bandwidth  will  cause  a 5.7  psec  incident  pulse  to  be  compressed  by 
a factor  of  T^/T  = 0.5  in  width.  Clearly,  larger  values  of  X2  are  required 


w. 
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to  extend  this  method  to  the  compression  of  CO2  laser  pulses  of  the 
order  of  0.1  nsec  in  width.  Calculations  show  that  significantly 
larger  values  of  cannot  be  realized  in  practice  by  increasing  the 
average  carrier  frequency  above  925  cm"\  because  the  attentuation 
length  decreases  rapidly  as  approaches  the  TO  phonon  frequency  at 
1072  cm"\ 

It  is  reasonable  to  inquire  whether  the  quadratic  dispersion 
parameter  might  be  significantly  increased  by  cooling  the  crystal  in 
order  to  decrease  anharmonic  damping  of  the  TO  phonon  and  thereby  obtain 
a larger  attenuation  length.  However,  at  lower  temperatures,  with 
correspondingly  smaller  damping  rates  the  curves  for  Re[k(w)] 
become  more  nearly  linear  at  frequencies  for  which  the  attenuation  length 
may  be  suitably  large.  Calculations  indicate  lRe(i2^)|^^^  can  be  in- 
creased by  only  a factor  of  about  1.3  at  frequencies  from  925  cm'^  to 
1060  cm"^  upon  cooling  the  crystal  from  300°K  to  80°K. 

Additionally,  it  should  be  noted  that  the  magnitude  of  varies 
significantly  over  the  large  THz  chirp  bandwidths  considered  in  the 
present  example.  Further  analysis,  to  be  published  separately,  reveals 
that  this  higher  order  variation  of  k(oj)  with  frequency  over  the  incident 
pulse  bandwidth  gives  an  asymmetrically  shaped  compressed  pulse  of  slightly 
greater  width  than  calculated  for  the  above  examples. 

SURFACE  PHONON -POLARITON  MODES  ON  BeO: 

Experimental  measurements  have  shown  that  the  attenuation  length  for 
a surface  polariton  may  be  several  times  greater  than  its  wavelength. 
Therefore,  surface  propagation  modes  may  also  be  considered  for  dispersive 
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optical  pulse  delay  lines.  The  theory  for  scattering  of  surface  polaritons 

by  surface  roughness  has  shown  that  the  quality  of  surface  required  to 

reduce  scattering  to  acceptable  limits  is  well  within  the  capabilities  of 

42 

optics  polishing  technology.  Theoretical  work  on  coupling  an  optical 

pulse  to  the  surface  mode  by  using  grating  structures  has  yielded  criteria 

for  launching  a frequency-chirped  surface  polariton  pulse  without  distorting 

43  44 

Its  shape  or  time-dependent  phase.  ’ 

Single-Interface  Structure:  A surface  optical  wave  may  propagate 

nonradiatively  along  the  interface  between  two  materials  having  dielectric 

constants  £p  and  provided  the  real  part  of  the  dielectric  constant  of 

one  medium,  Re(e^),  is  negative  and  greater  in  magnitude  than  the  positive 

5-7 

dielectric  constant  Cp  of  the  transparent  medium.  This  condition 
may  be  satisfied  for  frequencies  just  above  the  TO  phonon  frequency  of 
an  infrared-active  crystal.  The  dispersion  relation  for  a surface  polariton 
IS  of  the  same  form  as  given  by  Eq.  (10),  if  k(oj)  is  understood  to  be  the 
wave  vector  parallel  to  the  surface  and  the  "effective  dielectric  constant" 
£(oj)  for  the  surface  mode  is  defined  as: 
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axis  of  the  BeO  crystal  is  assumed  to  be  perpendicular  to  the  surface  wave 
vector  k and  the  electric  vector  E,  so  that  parameters  for  the  ordinary 

'\y 

ray  in  BeO  are  appropriate  for  calculating  £^(0)).  It  is  apparent  from 

the  slopes  of  the  solid  curves  that  group-velocity  dispersion  may  be  made 

large  over  any  desired  range  of  frequencies  just  above  ojjq  by  appropriate 

choice  of  the  transparent  medium  of  dielectric  constant  Cp.  For  example, 

if  £p=2.59,  the  surface  mode  exhibits  large  group-velocity  dispersion 

for  a range  of  frequencies  coincident  with  the  strong  (00°1 )-(10°0)  band 

P branch  emission  lines  from  a CO2  laser. 

Analytic  expressions  derived  from  Eqs.  (10)  and  (12)  for  the 
2 2 

derivatives  dk/dw  and  d k/doj  were  used  to  calculate  the  expansion  parameters 

and  T2  for  k(aj)  of  the  surface  mode.  The  frequency-dependent  dielectric 

constant  Ej(cj)  for  the  ordinary  ray  in  BeO  was  calculated  from  Eq.  (11), 

using  e(®)=2.95,  S=3.94,  and  oipQ=725.  cm"\^^  The  solid  curve  of  Fig.  5 

illustrates  the  dependence  of  T2  upon  average  carrier  frequency  for  the 

surface  phonon-pol ari ton  on  a Be0:air  surface,  with  damping  neglected. 

As  expected,  diverges  as  the  average  carrier  frequency  asymptotically 

approaches  the  maximum  frequency  oj^p  of  the  surface  mode. 

When  absorption  in  the  BeO  substrate  is  included,  the  solid  dispersion 

curve  in  Fig.  4 is  replaced  by  the  dashed  curve  for  Re[k(o))]  and  the  dotted 

curve  for  Im[k(o))].  The  damping  rate  y=7.5  cm'^  was  used  in  Eq.  (11)  to 
45 

calculate  The  curvatures  of  the  broken  curves  in  Fig.  4 indicate 

2 

that  real  and  imaginary  components  of  T2  display  maxima  at  frequencies 
close  to  the  limiting  surface  frequency  oj^p.  Therefore,  as  for  bulk 
phonon-pol ari tons , the  effects  of  absorption  in  the  BeO  substrate  are  to 
limit  both  the  maximum  attainable  value  of  T2  and  the  attenuation  length. 
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2 

The  dependence  of  components  of  i2  upon  average  carrier  frequency  is 
illustrated  in  Fig.  5 for  a BeO:air  structure  with  x=l  cm.  The  dashed 

curve  displays  maxima  for  positive  values  of  Re{x2  ) at  1017  cm"  and 

_ 1 2-1 
1028  cm  frequency,  and  a peak  negative  value  for  Re  (t2  ) at  1024  cm 

2 

frequency.  Similarly,  the  dotted  curve  for  Im(T2  ) displays  positive 
maxima  at  1022  cm  ^ and  1032  cm'^  and  a peak  negative  value  at  1026  cm~^ 
frequency. 

As  an  example,  consider  a pulse  having  an  average  carrier  frequency 

of  980  cm~\  corresponding  to  the  (00°1 )-(10°0)  band  R(26)  emission  line 

of  a CO2  laser,  for  which  the  attenuation  length  of  a BeOrair  surface  mode 

is  20  ym, (Re(x2^) 1 2x10'^^  for  x = 1cm,  and  Im(x2^)=0. 23Re(x2^) . 

Propagation  a distance  of  0.2  mm  then  gives  a quadratic  dispersion  parameter 
2 1/2  «1 2 

lRe(x2  )|  ' =2.8x10"  sec,  with  a reduction  in  peak  intensity  by  a factor 
exp[-2Im(k^)x]  = 4x10"^  due  to  absorption.  As  was  the  case  for  the  bulk 
phonon-polariton  in  quartz,  a quadratic  dispersion  parameter  of  this 
magnitude  requires  that  the  initial  chirp  bandwidth  be  a few  THz  and  the 
incident  pulse  width  of  the  order  of  a few  picoseconds.  If,  for  example, 
the  frequency  chirp  covers  a range  of  -0.9  THz,  corresponding  to  the  entire 
(00°1 ) -( 10°0)  band  R branch  of  a highly  pressure  broadened  CO2 
laser,  then  an  11  psec  pulse  may  be  optimally  compressed  to  0.7  psec  in 
width,  with  gain  in  peak  intensity  by  a factor  of  T/T^=16.  A 22  psec 
pulse  will  be  compressed  by  a factor  T^/T=0.5,  using  the  same  initial  chirp 
bandwidth.  As  previously  noted,  additional  pulse  shaping  is  anticipated  due 
to  variation  of  T2  with  frequency  over  the  large  chirp  bandwidth.  However, 
it  is  clear  that  pulse  compression  is  limited  by  the  magnitudes  of  T2  and 
Im(kQ)  to  incident  pulses  of  the  order  of  20  psec  in  width. 
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Layered  Structure:  Another  configuration  of  interest  is  a layered 

structure  formed  by  depositing  a thin  film  of  transparent  dielectric  onto 
a surface-active  material  such  as  BeO.  The  surface  polariton  in  this 
layered  structure  has  a dispersion  relation  of  the  form  given  by  Eq.  (10), 
but  for  this  case  the  "effective  dielectric  constant  e(u)"  for  the  surface 
mode  is  given  by  the  root  of  the  following  transcendental  equation:^^ 


(13) 


where  Cp,  and  e^(oi)  are  the  dielectric  constants  for  the  semi-infinite 
transparent  medium  (air),  the  thin  dielectric  overlayer  film  of  thickness  l, 
and  the  surface-acti medium  (BeO),  respectively.  The  detailed  configura- 
tion and  expressions  for  the  surface  polariton  fields  are  outlined  in 
reference  46,  along  with  the  derivation  of  the  transcendental  equation. 

We  are  presently  interested  in  that  portion  of  the  dispersion  curve  for 
which  £p<e<e^,  for  which  the  electric  field  amplitude  varies  sinusoidally 
through  the  film  region  in  the  direction  perpendicular  to  the  surfaces. 

When  absorption  is  neglected,  the  solid  dispersion  curve  in  Fig.  4 
increases  to  a maximum  allowed  frequency,  then  bends  downward  for  larger 
values  of  surface  wave  vector  k.  At  frequencies  near  the  maximum  in  the 
dispersion  curve,  large  values  for  X2  occur  at  relatively  small  values  of 
wave  vector.  Therefore,  it  should  be  possible  to  couple  an  optical  beam 
to  the  surface  mode  at  a point  on  the  dispersion  curve  where  k(aj)  varies 
rapidly  with  frequency,  using  either  the  prism-coupling  method  or  a grating 


r 
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coupler  with  reasonably  large  period.  The  layered  structure  also  offers 
the  potential  advantage  that  both  the  dielectric  constant  and  the 
thickness  of  the  overlayer  film  are  experimental  variables  which  may 
be  used  to  optimize  the  design  of  a dispersive  pulse  delay  line. 

Figure  4 illustrates  the  dispersion  curve  for  a layered  structure 
consisting  of  a germanium  film  on  a BeO  crystal,  with  a germinium  film 
thickness  of  500  A.  The  c axis  of  the  Beo  crystal  is  assumed  parallel  to 
the  plane  of  the  interface  and  perpendicular  to  the  direction  of 
propagation,  so  that  parameters  for  the  ordinary  phonon-polariton  of 
BeO  are  used  to  calculate  £^(0)).  The  solid  curve  is  calculated  for 
the  limit  of  no  absorption  in  the  BeO  substrate.  Similar  calculations 
show  that  a thicker  Ge  film  lowers  the  maximum  allowable  frequency  and 
causes  the  dispersion  curve  to  bend  more  sharply  in  the  vicinity  of 
this  maximum.  The  solid  dispersion  curve  bends  downward  and  asymptot- 
ically approaches  from  above  the  limiting  surface  frequency  for  a 
single-interface  structure  of  BeO:Ge  at  very  large  values  of  wave  vector. 

The  dispersion  parameters  for  a surface  polariton  in  a layered 
structure  are  calculated  from  analytic  expressions  derived  from  Eqs.  (10) 
and  (13).  Formulation  of  expressions  for  de/du  and  d e/do  from  Eq.  (13) 
is  much  less  straightforward  than  for  the  previous  configurations,  so 
that  the  calculation  for  the  layered  structure  deserves  special  comment. 

It  should  be  noted  that  the  implicit  dependence  of  e upon  frequency 
through  £^(0)  must  be  included,  in  addition  to  the  explicit  dependence 
upon  o in  the  argument  of  tan(o£(£^-£)^'^^/c)  in  Eq.  (13).  Further  details 
of  the  computation  are  deferred  to  Appendix  A,  since  the  expressions  are 
lengthy  and  do  not  contribute  insight  to  the  present  discussion. 


11-18 


Variation  of  with  average  carrier  frequency  is  illustrated  in 
Fig.  6 for  Ge  overlayer  thicknesses  t=100  ^ andt=500  A and  x=l  cm. 

The  case  for  no  absorption  (y=0)  is  not  represented  in  Fig.  6 because, 
as  is  evident  from  the  solid  dispersion  curve  in  Fig.  4,  the  frequency 
dependence  of  k(u),  and  therefore  ofi^,  is  double-valued.  However, 
it  is  clear  from  Fig.  4 that  for  no  damping  i2  diverges  at  frequencies 
approaching  the  region  of  zero  slope  in  the  solid  dispersion  curve,  at 
943  cm"^  fort=500  A.  Positive  values  of  i2  occur  for  wave  vector  values 
below  this  point,  and  negative  values  of  T2  for  larger  values  of  k(oj). 

When  absorption  in  the  BeO  region  is  included.  Fig.  6 shows  that 

I Re(T2^) 1 for  the  100  A thick  film  displays  maxima  in  its  frequency 

2 -11 
dependence,  with  Re(t2  )<0  at  frequencies  between  995  cm"  and  1014  cm" 

2 2 

and  Re(t2  )>0  otherwise;  while  Im(T2  ) is  negative  for  frequencies  above 

999  cm'^  and  positive  below  999  cm"\  The  500  A thick  Ge  film  gives 

2 1 -1 
similar  structure,  with  Re(t2  )<0  between  940  cm"  and  983  cm  and 

2 -1 

Re(t2  )>0  above  945  cm"  frequency.  The  maximum  value  of  T2  attainable 

for  the  layered  structure  is  about  an  order  of  magnitude  smaller  than  for 

the  single-interface  surface  polariton  on  Be0:air  or  the  bulk  phonon- 

polariton  in  a -quartz,  because  of  the  influence  of  absorption. 

As  an  example,  consider  a pulse  having  an  average  carrier  frequency 

of  900  cm"\  corresponding  to  the  (01°l)-(03^0)  band  P(31)  emission  line  of  a 

CO2  laser,  propagating  as  a surface  polariton  in  a Be0:Ge  filmrair  structure 

with  a film  thickness  )l=100  A.  The  attenuation  length  for  this  frequency 

is  55ym  and,  for  x=l  cm,  | Re(t2^) | ^'^^=6x10"^^  sec,  with  |Im(T2^)l^^^  = 

2 1/2 

0.07  |Re(Tp  )|  ' . A propagation  distance  of  0.63  mm  then  gives 
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2 1/2  - 1 2 

|Re(x2  ) ' =1.5x10"  sec,  along  with  an  attenuation  factor 

_ 5 

exp[-2Im(k^)x]=10  due  to  absorption.  Therefore,  for  the  layered 
structure  is  comparable  to  values  obtained  with  similar  attenuation 
for  a single  BeOrair  surface,  and  pulse  compression  is  limited  to  pulses 
of  a few  picoseconds  initial  width  using  THz  chirp  bandwidths. 

SUMMARY  AND  CONCLUSIONS 

An  analysis  has  been  presented  within  the  framework  of  linear 

dispersion  theory  for  the  propagation  and  temporal  compression  of 

linearly  chirped  pulses  in  dispersive,  absorbing  media.  The  relative 

delay  among  frequency  groups  in  the  incident  pulse  is  characterized 

in  terms  of  the  frequency  dependence  of  wave  vector  k(w)  for  each  type 

of  propagation  mode  considered.  Specific  calculations  are  presented 

for  CO2  laser  pulses  which  propagate  as  bulk  phonon-polaritons  in 

quartz,  as  surface  phonon-polaritons  on  a BeO:air  surface,  and  as 

surface  phonon-polaritons  in  a BeO:Ge  filmrair  layered  structure. 

For  attenuation  losses  of  50  dB  or  smaller  significant  pulse 

compression  may  be  demonstrated  for  optical  pulses  20  psec  or  less 

in  width  having  applied  chirp  bandwidths  in  the  THz  range.  Larger 

group  dispersion  and  smaller  absorption  loss  are  required  for  this 

method  to  become  competitive  with  other  types  of  optical  pulse 

15  19-22 

compression  filters,  such  as  the  grating  pair  ’ or  a Gires- 

23  24 

Tournois  interferometer,  ’ or  to  extend  this  method  to  the  com- 
pression of  nanosecond  CO2  TEA-laser  pulses  using  experimentally 
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realizable  chirp  bandwidths.  Appropriate  materials  must  have  larger 
TO  Dhonon  oscillator  strengths  and  smaller  anharmonic  damping  rates 
at  the  appropriate  frequencies.  Further,  the  limitations  imposed  by 
laser-induced  damage  upon  peak  and  average  laser  power  should  be  con- 
sidered when  assessing  particular  applications.  However,  if  a suitably 
dispersive  medium  can  be  identified,  the  simplicity  and  compactness  of 
this  type  of  pulse  compression  filter  may  become  important  design 
considerations . 
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APPENDIX  A 

COMPUTATIONAL  METHOD  FOR  LAYERED  STRUCTURES 

Analytic  expressions  for  the  propagation  constants  for  bulk  phonon- 
polaritons  and  for  surface  phonon-polaritons  on  a single-interface 
structure  may  be  derived  in  a straightforward  manner  from  the  dispersion 
relation  for  k{aj).  However,  for  the  layered  structure  the  propagation 
constants  must  satisfy  the  transcendental  equation  (13),  so  that  the 
development  of  analytic  expressions  for  is  much  less  straightforward. 
Therefore,  the  method  used  to  calculate  k(w)  and  its  frequency  derivatives 
for  layered  structures  is  briefly  outlined  here. 

The  transcendental  equation  (13)  is  of  the  general  form 

F[c(lo,  £^(oi))]  = 0 (A-1) 

The  wave  vector  k(ijj)  and  its  frequency  derivatives  are  calculated  using 
the  complex  root  c(tj),  for  which  both  real  and  imaginary  parts  of  F(e) 
equal  zero.  The  roots  e(u)  are  conveniently  found  by  mapping  the  root 
locus  F(c)=0  in  the  complex  e plane.  The  root  e(u))  is  calculated  first 
at  a frequency  S^Q+fioj  just  above  the  TO  phonon  frequency,  where  Re(e^)  is 
large  in  magnitude  and  c((jJjq+6w)  is  very  nearly  equal  to  the  value 
calculated  directly  from  Eq.  (12)  without  the  thin  overlayer  film. 

Equation  (12)  provides  a first  approximation  for  e(ajjQ+6w),  valid  for 
very  thin  dielectric  overlayers.  Based  upon  this  initial  estimate,  the 
exact  root  at  wjq+6oj  is  calculated  by  a Newton-Raphson  iteration  in  complex 
variables,  using  analytic  expressions  formed  for  F(c)  and  dF/de.  The 
frequency  is  then  incremented  to  u)yQ+26w,  usually  by  Su)  = 5 cm'^  or  less. 
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and  e(u))  is  calculated  for  the  new  frequency  using  as  an 

initial  value  for  the  Newton-Raphson  iteration.  In  this  manner,  the 
root  locus  of  solutions  e(u))  is  traced  out  step-wise  over  the  entire 
range  of  surface  polariton  frequencies,  and  each  solution  rapidly 
converges  in  a small  number  of  iterations. 

Calculation  of  T2  requires  the  derivation  of  expressions  for  the 
2 2 

derivatives  dc/doj  and  d e/cku  . The  transcendental  function  F(e)  is 
a function  of  frequency  through  e((i))  which,  in  turn,  varies  with 
and  also  explicitly  with  uj  in  the  argument  of  tan(aj)!,(e^-e)  ' /c) . If 
the  partial  derivative  with  respect  to  e is  denoted  by  a prime  (') 
notation  and  the  partial  derivative  with  respect  to  frequency  by  a dot 
(■)  notation,  then  the  frequency  derivative  of  F is  written; 


(A-2) 


which  yields  a direct  analytic  expression  for  de/du.  Differentiation  a 
second  time  gives: 


^ = F + 2 F'  ^ + F"  f + F‘  ^ = 0 

.4  dca  { QUi  I .4 

cjd  y / du) 


(A-3) 


2 2 

which  similarly  yields  a direct  analytic  expression  for  d e/doj  . The  par- 
tial derivatives  F'  and  F"  are  formed  from  the  explicit  appearance  of  e in 
Eq.  (A-1),  while  derivatives  F and  F must  include  the  frequency  dependence 
of  as  well  as  an  explicit  dependence  on  frequency. 
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Having  now  formed  analytic  expressions  for  e(w),  de/doj  and  d e/dco  , 
the  propagation  constant  k(oj)  and  its  frequency  derivatives  are 


? 1 


i 


f 1 


evaluated  from  Eq.  (10)  to  give: 
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II-l.  Pulse-width  compression  ratio  T^/T  versus  normalized  quadratic  disper- 
sion parameter  x-  for  various  ratios  of  6ui  /6tii  . Dashed  curves  ( ) 

show  only  pulse  broadening  for  while  solid  curves  illustrate 

pulse  compression  when  Optimum  pulse  compression,  illustrated 

2 

by  the  dotted  curve  (••••) i occurs  when  5u^/5cj^=-(t2/T)  . 

II-2.  Dispersion  curves  for  ordinary  (^c)  bulk  phonon-polaritons  in  a -quartz. 
Solid  curve  gives  dispersion  in  the  limit  of  no  anharmonic  damping 

Y|^=0;  and  broken  curves  give  Real  [k(u))/2Ti]  ( ) and  Im[k(aj)/2iT]  (••••) 

at  300°K. 


II-3.  Variation  with  average  carrier  frequency  of  the  quadratic  dispersion 

parameter  T2  for  the  ordinary  bulk  phonon-polariton  of  a -quartz, 

normalized  to  a propagation  distance  x.l  cm.  Solid  curve  illustrates 

the  limit  for  no  damping,  while  the  components  when  absorption  at 

300°K  is  included  are  shown  for  | Re(t2^)  I by  the  dashed  curve  ( ) 

2 1/2 

and  for  |lm(x2  )|  ' by  the  dotted  curve  (••••)■ 

II-4.  Dispersion  curves  for  ordinary  (Elc)  surface  phonon-polaritons  on  single- 

interface  structures  of  BeO;air  (e  =1.)  and  BeOrCsI  (e  =2.59),  and  on  a 

P P 

layered  structure  of  BeOrGe  film:air  with  Ge  film  thickness  £=500  A. 

Solid  curves  are  for  the  limit  of  no  anharmonic  damping.  With  absorp- 
tion in  BeO  included,  components  of  complex  k(uj)  are  indicated  by  the 

dashed  curves  ( ) for  Re[k(oj)]  and  the  dotted  curves  (••••)  for 

Im[k(oj)]. 

II-5.  Variation  with  carrier  frequency  of  the  quadratic  dispersion  parameter 
T2  for  the  ordinary  surface  phonon-polariton  on  a BeOiair  interface. 


nonnalized  to  a propagation  distance  x=l  cm.  Solid  curve  illustrates 

the  limit  for  no  damping,  while  the  components  when  absorption  at 

300°K  is  included  are  shown  for  |Re(T2  )|  ' by  the  dashed  curve  ( ) 

2 1/2 

and  for  1 Im(Tp  ) 1 by  the  dotted  curve  (••••)• 

II-6.  Variation  with  carrier  frequency  of  the  quadratic  dispersion  parameter 
for  the  ordinary  surface  phonon-polari ton  in  layered  structures  of 
8eO;Ge  filmrair,  with  Ge  film  thicknesses  Jl=100  A and  £.=500  A. 

2 1/2 

Absorption  in  the  BeO  region  is  included,  with  components  iReCip  )|  ' 

2 1/2 

indicated  by  the  dashed  ( ) curves  and  | Im(T2  )1  ' by  the  dotted 

(....)  curves. 
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Ill,  Chirped  Pulse  Compression  In  Strongly  Dispersive  Media* 


J.  D.  McMullen 
Electronics  Research  Center 
Rockwell  International  Corporation 
Anaheim,  California  92803 


ABSTRACT 


An  analysis  is  given  for  the  temporal  compression  of  linearly 
chirped  pulses  in  strongly  dispersive  media.  Terms  cubic  in  fre- 
quency are  included  in  the  frequency  dependent  wave  vector  k(u))  of 
the  dispersive  medium.  Analytical  expressions  are  developed  for  the 
compressed  pulse  envelope,  and  numerical  examples  given  for  the 
influence  of  the  cubic  frequency  contribution.  The  strong  dispersion 
term  is  found  to  give  asymmetric  broadening  of  the  compressed  pulse 
envelope. 


♦This  research  was  supported  by  the  U.  S.  Army  Research  Office, 
under  Contract  No.  DAHC04-74-C-0024. 
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Introduction 

A pulse  with  a linearly  chirped  carrier  frequency  may  be  temporally 

1-13 

compressed  upon  propagating  through  a suitably  dispersive  medium. 

Previous  analyses  have  been  given  for  chirped  pulses  propagating  in  dis- 

1-5 

persive  but  nonabsorbing  media,  and  for  chirped  pulses  in  absorbing, 
moderately  dispersive  media. ^ These  treatments  are  valid  when  the  dis- 
persion is  sufficiently  weak  that  the  group  velocities  change  only  linearly 
with  frequency  across  the  bandwidth  of  the  incident  pulse. 

A recursion  solution  has  been  reported  for  chirped  pulses  propagating 

through  strongly  dispersive  but  nonabsorbing  media, ^ and  for  non-chirped 

7-9 

pulses  in  strongly  dispersive  and  absorbing  media.  Numerical  calcu- 
lations of  pulse  shapes  have  been  reported  for  non-chirped  pulses  in 

dispersive,  absorbing  media,  requiring  evaluation  by  computer  of  the  in- 

13 

verse  Fourier  transform  for  the  pulse  shape.  For  situations  where  the 
strong  dispersion  results  from  an  average  carrier  frequency  lying  close 
to  an  absorption  frequency  of  the  medium,  frequency-dependent  absorption 
can  play  an  important  role  in  determining  the  shape  and  final  width  of 
the  compressed  pulse. 

In  this  paper,  closed-form  analytic  solutions  are  presented  for  a 
linearly  chirped  pulse  propagating  in  an  absorbing,  strongly  dispersive 
medium,  where  terms  through  cubic  in  frequency  are  retained  in  the  fre- 
quency dependence  of  the  wave  vector  k(o)).  Numerical  examples  are  given 
to  illustrate  the  contribution  of  this  cubic  dispersion  term  to  the  shape 
of  the  output  compressed  pulse  envelope,  and  the  limitations  imposed  by  it 
upon  the  attainable  compressed  pulse  width. 
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Theory 

Consider  a linearly  chirped  pulse  having  average  carrier  frequency 
(rad/sec)  and  a Gaussian  envelope  of  the  form: 


specified  at  x=0  just  inside  the  boundary  of  a frequency-dispersive 

-1  2 

medium.  The  full  width  at  e’  of  peak  intensity  |E^|  is  T,  and  the  range 
of  frequencies  swept  over  during  the  full  width  is  (rad/sec).  The 
Fourier  transform  of  E(0,t)  is  given  by: 


f 2 2 1 ^ 

6w  = ('5Wq)  + 


(3) 


-1  i~  2 

is  the  full  width  at  e of  peak  for  the  power  spectrum  |E(0,w)|  , and 

= 4/T  is  the  bandwidth,  similarly  defined,  for  the  power  spectrum  of 

the  unmodulated  pulse  envelope.  The  pulse  is  assumed  to  propagate  in  the 

+x  direction  in  the  dispersive  medium,  which  is  characterized  by  complex 

wave  vector  k(w).  If  the  first  term  of  E(0,oj)  in  Eq.  (2)  peaked  in  ampli- 

. ( + ) 

tude  about  w=cj  is  denoted  as  E{0,.)  and  the  second  term  peaked  at  co=-u  as 
E(0,uj),  the  output  pulse  waveform  ;t  position  x is: 
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E(x.t) 


'/j*“  ( 


^ (+) 
E(0,w) 


ik('ji))x 


^ (-) 
E(O.u)) 


■1k(oj)x 


-iwt 


(4) 


which  ensures  that  E(x-+®,t)  ->  0 in  the  presence  of  absorption  in  the  medium. 
The  dispersive  medium  is  characterized  by  a complex,  frequency-dependent 
wave  vector  k(cj)  which  is  assumed  to  be  slowly  varying  with  w over  the  band- 
width of  the  incident  pulse.  The  wave  vector  k{oj)  is  expressed  by  the 
Taylor  series 

1 / ^ T 

k(oj)x  = k^x  + Ti(w-u)^)  5"  ( 2~)  +....  (5-a) 

about  aj=+u^,  and  (5-a) 

k(w)x  = k^x  - *r-j(a)+cu^)  + ^ ~ +....  (5-b) 

about  a)=-cjp,  where  k^  5 expansion  coefficients  are  defined 

as 


_ dk 
^1  = ^ 


= 4 


d^k 


and 


3 = i d-"k 


(5-c) 


(5-d) 


(5-e) 


The  purpose  of  this  form  of  definition  for  and  is  to  obtain  simple 
parameters  of  the  dispersive  medium,  having  units  of  time,  which  can  be 
compared  to  the  incident  pulse  width  and  bandwidth  for  the  discussions  of 
pulse  compression  to  follow. 


'■"♦v  “ ^ . ’V  > i** 
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The  output  pulse  waveform  at  position  x is  evaluated  to  be: 


where  for  convenience  we  define 


2 


which  contains  information  about  the  initial  pulse  width,  the  applied 

chirp,  and  the  quadratic  dispersion  parameter  xo,  Ai(  ) is  the  well-knov/n 
14  ^ / 

Airy  integral.  Details  of  the  derivation  of  Eq.  (6)  are  given  in  Appendix 

A. 

A case  of  particular  interest  occurs  when  is  sufficiently  small 
that  the  cubic  term  in  the  expansion  for  k(to)  contributes  only  sMghtly 
to  shaping  of  the  output  pulse  envelope,  or  | -t)] j «| a^| , for 
which  E(x,t)  reduces  to  the  following  form: 


(8) 


r 


i 
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where  T is  defined  by: 

X 


(9) 


and  the  contribution  of  the  cubic  frequency  term  in  k{uj)  is  contained 
entirely  within  the  shape  factor  S(x,t),  given  by: 


S(x,t)  = exp 


3 

(2a)^  . 


(10) 


2 

In  the  limit  that  ^2=0,  the  shape  factor  S becomes  unity  and  |E(x,t)| 
reduces  to  a Gaussian  pulse  with  full  width  at  e"^  of  peak  intensity. 
For  this  case  Eq.  8 shows  the  peak  pulse  intensity  is  reduced  or  increased 
by  a factor  (T/T^)exp(-2Im[kQx])  as  the  output  pulse  becomes  either 
dispersion  broadened  or  compressed  in  width.  When  the  applied  chirp  is 
an  optimum  value  given  by 


then  the  output  pulse  E(x,t)  is  compressed  to  its  narrowest  possible  width. 

Figure  1 illustrates  the  compressed  pulse  waveform  for  a chirped  Gaussian 
pulse  having  initial  width  T = 1 nsec  and  an  applied  chirp  5a)^=-1006(jjQ  after  the 
pulse  has  propagated  through  a dispersive  medium  for  which  lm(k^x)=0  and 
T2=0.1  nsec.  Both  Eqs.  (6)  and  (8)  were  used  to  calculate  pulse  envelopes. 

The  approximate  solution  of  Eq.  (8)  greatly  simplifies  the  calculations  for 
ranges  of  parameters  such  that  the  argument  of  the  Airy  function  is  very 
large.  The  solid  curve  for  the  case  i2=0  shows  the  compressed  pulse  has  a 
Gaussian  profile  of  full  width  T^=10  psec  at  e"^  of  peak  intensity.  The 


4.. 


J 
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broken  curves  illustrate  the  influence  of  the  cubic  term  in  k(io)  upon  the 
output  pulse  shape  for  successively  larger  values  of  with  the  other 
parameters  held  constant.  For  larger  values  of  the  peak  of  the  compressed 
pulse  is  diminished  in  height  and  shifted  to  larger  delay  times,  and  the 
pulse  shape  becomes  asymmetrically  broadened  about  the  peak.  Interference 
structure  occurs  in  the  trailing  edge  of  each  pulse  for  T2=0  in  accordance 
with  the  oscillatory  nature  of  the  Airy  function  for  negative  real  argument. 

It  is  apparent  from  Figure  1 that  significant  pulse  shaping  effects 
occur  for  as  small  as  3 psec  in  the  above  example.  Upon  replacing  the 
factor  (w-Wp)  in  Eg.  (5-a)  by  the  bandwidth  of  the  incident  pulse,  it 
is  noted  that  the  magnitude  of  the  third  order  term  in  the  series  expansion 
for  k((ij)  is  only  about  0.3  percent  of  the  magnitude  of  the  second  order 
term.  Therefore,  the  influence  of  upon  the  shape  and  height  of  the 
compressed  pulse  becomes  quite  significant  for  smaller  values  of  than 
might  be  anticipated  upon  examination  of  Eg.  (5-a)  alone. 

Conclusion 

In  conclusion,  it  has  been  shown  that  a closed-form  analytic  solution 
may  be  obtained  for  a linearly  chirped  pulse  which  becomes  temporally 
compressed  upon  propagating  through  a strongly  dispersive  medium,  where 
terms  cubic  in  frequency  are  retained  in  the  expansion  of  k(w)  about  the 
average  carrier  frequency.  The  cubic  frequency  dependence  is  shown  to 
further  delay  the  peak  of  the  compressed  pulse  and  to  cause  the  pulse  to 
become  asymmetrically  broadened  about  that  peak.  The  small  value  of 
required  to  influence  the  compressed  pulse  shape  suggests  may  easily  be 
significant  in  pulse  compression  experiments  for  which  the  group  delay 
changes  rapidly  with  frequency  over  the  bandwidth  of  the  incident  pulse. 
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Appendix  A 

Analysis  of  the  output  pulse  shape  from  Eg.  (4)  requires  the 
evaluation  of  inverse  Fourier  transform  integrals  of  the  form: 


” r T ^ 

I = f dy  expji(T^-t)y  - ay^  + i 

where  y=<i)-a)^  and  a is  defined  by  Eq. 
Upon  introducing  a change  in  variable 

« = y ^ i f-3 
3 


The  integral  becomes: 
I = exp 


(?[ 


, . 2 ^ 
3^3 
^3 


)/■ 


dE  exp 


A-1 


A-2 


A-3 


14 

The  integral  in  Eq.  (A-3)  is  recognizable  as  the  Airy  integral,  leading 
to  the  analytic  expression  in  Eq.  (6)  for  the  pulse  shape  at  position  x. 

When  Tj  is  very  small  and  the  argument  of  the  Airy  function  becomes 
very  large  in  magnitude,  it  is  useful  to  approximate  Ai  by  the  following 
1 imit: 


Lim  Ai(z) 

|z|->« 


where 


” T 

1 + (^1  -t) 

^ d ^ 


A-4 


A-5 


III-9 


Series  expansion  of  must  be  carried  out  to  include  terms  quadratic  in 
because  of  the  exponential  factor  proceeding  the  integral  in  En.  (A-3), 


so  that 


1 /4 

while  the  expansion  for  z in  the  denominator  of  Eq.  A-4  may  be  truncated 

3 

after  the  first  order  term  in  as  follows: 


.1/4 


J/2 


1+1  — 2 
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The  resulting  expression  for  the  output  pulse  shape  in  Eq.  (8)  greatly 
reduces  the  computational  effort  for  E(x,t)  for  ranges  of  parameters  over 
which  the  argument  of  the  Airy  function  is  large  in  magnitude. 
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Figure  Captions 

Figure  III-l  Waveform  of  compressed  chirped  pulse  intensity  with  T-1  nsec, 

=-1006w  , and  t-,=0.1  nsec,  for  various  values  of  group 
mot 

dispersion  1^=0  ( ),  T2=3  psec  ( ),  T2=4  psec  ( ), 

13=5  psec  { ) and  t2=10  psec  ( ). 
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IV.  THEORY  OF  SURFACE  POLARITON  GENERATION  THROUGH  GRATING  COUPLERS* 

D.  L.  Mills 
Department  of  Physics 
University  of  California 
Irvine,  California  92664 

ABSTRACT 

We  discuss  the  theory  of  surface  polariton  generation  by  grating  couplers. 
The  emphasis  is  on  an  analysis  of  the  manner  in  which  the  grating  may  modify  the 
frequency  spectrum  of  a shaped  incident  pulse.  If  L is  the  length  of  the  portion 
of  the  grating  illuminated  by  the  incident  beam,  then  we  find  the  frequency  spec- 
trum of  the  incident  pulse  can  be  transmitted  to  the  surface  polariton  undistor- 
ted if  the  criterion  L < V^/o  is  satisfied,  where  is  the  group  velocity  of  the 
surface  polariton,  and  a a measure  of  the  spread  of  frequencies  in  the  incident 
beam  used  to  excite  the  surface  polaritons.  If  this  criterion  is  not  satisfied, 
the  frequency  spectrum  of  the  transmitted  pulse  will  be  narrower  than  that  in  the 
incident  pulse.  The  treatment  here  provides  expressions  for  the  absolute  ampli- 
tude and  shape  of  the  transmitted  signal,  within  i , framework  of  a method  that 
assumes  the  height  of  the  grating  is  small. 
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A . Introductory  Remarks 

There  is  considerable  current  interest  in  the  study  of  surface  polaritons. 

These  are  electromagnetic  surface  waves^^^  that  under  a variety  of  circumstances  may 
propagate  along  an  interface  between  two  media,  with  electromagnetic  fields  confined 
to  the  near  vicinity  of  the  interface.  For  the  simple  case  of  the  interface  between 
vacuum  and  an  isotropic  dielectric  with  frequency  dependent  dielectric  constant  e(w), 
these  waves  may  propagate  whenever  Re(c(u))  < 0.  One  may  launch  them  through  use  of 
a prism  coupling  device  or  a grating^ ^ ^ and  their  propagation  length  may  be  in 
the  range  of  1 mm  to  1 cm  at  10. 6u,  for  selected  substrates. 

Since  the  propagation  length  of  these  waves  can  be  reasonably  long  in  the  infra- 
red, in  principle  one  may  contemplate  incorporating  them  into  devices.  There  can  be 
substantial  curvature  in  the  surface  polariton  dispersion  relation,  and  as  a conse- 
quence one  may  consider  the  use  of  surface  polaritons  for  signal  processing  opera- 

(31 

tions  such  as  compression  of  a chirped  input  pulse;  'Experimental  studies  of  this 
possibility  are  underway,  but  they  are  currently  in  an  early  stage. 

In  such  an  operation,  one  envisions  a chirped  pulse  incident  on  a coupler  which 
launches  a surface  polariton  pulse  that  propagates  on  a substrate.  One  is  then  led 
to  inquire  into  the  question  of  whether  the  surface  polariton  pulse  is  composed  of  a 
frequency  spectrum  that  accurately  reflects  that  in  the  input  pulse.  Since  the 
response  function  of  the  grating  coupler  is  sharply  peaked  about  a particular  fre- 
quency (that  of  a surface  polariton)  when  excited  by  radiation  incident  at  a fixed 
angle,  the  question  seems  particularly  important  to  explore  in  this  case. 

The  purpose  of  this  paper  is  to  present  a study  of  the  shape  of  the  surface 
polariton  pulse  launched  through  use  of  a grating  coupler  excited  by  an  incident 
wave  with  a spectrum  of  frequencies  contained  in  it.  We  consider  the  specific 
case  of  a chirped  input  pulse,  and  find  expressions  for  the  absolute  amplitude  and 
shape  of  the  output  surface  polariton  pulse,  when  the  height  of  the  grating 


L 


Sit- 


i 


i 


i 


segment  illuminated  by  the  incident  beam,  is  the  group  velocity  of  the  surface 
polariton  launched  from  the  coupler,  and  o the  width  in  frequency  of  the  input  beam 
(in  radians/sec),  then  provided 


L < Vg/a  , 


(A-1) 


the  frequency  composition  of  the  surface  polariton  generated  can  accurately  reflect 
that  of  the  input  pulse.  When  one  has  instead  L > Vg/o,  the  frequency  spectrum  of 
the  output  pulse  is  narrower  than  of  the  input  pulse.  This  is  most  unfavorable, 
since  the  coupler  then  acts  to  expand  the  surface  polariton  pulse,  thus  making  the 
task  of  subsequently  compressing  the  pulse  more  difficult. 

In  essence,  the  criterion  quoted  above  requires  one  to  focus  the  incident  beam 
down  to  a size  small  compared  to  Vg/o,  if  one  wishes  to  pass  the  frequency  spectrum 
in  the  incoming  pulse  faithfully  on  to  the  surface  polariton  pulse  generated  by  the 
grating.  Insertion  of  numbers  into  the  criterion  indicates  that  the  focusing 
requirement  places  no  serious  limitation  on  the  use  of  gratings,  although  for  pulses 
that  approach  picosecond  duration,  one  may  well  need  to  worry  about  this  question. 

The  organization  of  this  paper  is  as  follows.  In  Section  B , we  derive  an 
expression  for  the  amplitude  of  the  surface  polariton  produced  by  a linearly  chirped 
pulse  incident  on  a grating.  This  derivation  leans  heavily  on  results  obtained  in 
an  earlier  paper  by  the  present  author.  The  result  obtained  in  Section  B is  of 
reasonable  generality,  and  after  some  manipulation  we  cast  it  into  a form  that  is 
useful  for  numerical  computation.  In  Section  C , we  work  out  the  form  of  the  sur- 
face polariton  pulse  for  several  special  limiting  situations.  These  results, 
obtained  in  analytic  form,  allow  the  principal  parameters  of  the  problem  to  be 
isolated,  and  their  role  evaluated. 
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B . Derivation  of  an  Expression  for  the  Shape  of  the  Output  Pulse 

The  geometry  that  forms  the  basis  of  the  present  discussion  is  illustrated  in 
Figure  (1).  An  incident  beam  of  finite  width  is  presumed  to  strike  a diffraction 
grating  ruled  on  a substrate  with  dielectric  constant  e(u).  The  angle  of  incidence 
is  Op,  and  to  generate  surface  polaritons  the  incident  radiation  is  p-polarized. 

The  plane  of  incidence  of  the  incident  radiation  is  perpendicular  to  the  grooves  in 
the  grating. 

Consider  the  incident  wave  to  be  monochromatic  with  frequency  w for  the  moment 
and  let  d be  the  grating  spacing.  The  wave  vector  component  of  the  incident  radia- 
tion parallel  to  the  surface  is  usino^  /c,  with  c the  velocity  of  light. 

Near  the  surface  of  the  grating,  the  electromagnetic  field  of  the  incident  wave  is 
perturbed  from  the  plane  wave  form,  and  acquires  ripples  with  wave  vector  parallel 
to  the  surface 

^ sin  Op  + ^ n ; n = 0,  1,  2 ...  . (B-1) 


Derivation  of  the  detailed  form  of  the  field  near  the  grating  has  been  presented  in 
recent  works  by  Elson^^^  and  by  the  present  author. (While  the  main  emphasis  in 
Reference  (6)  differs  from  the  question  explored  here,  a prescription  is  given 
there  for  constructing  the  fields  near  the  grating,  and  it  is  this  prescription  we 
have  used  here.) 

Now  let  Wg(k„)  be  the  dispersion  relation  of  the  surface  polariton.  For  the 
semi-infinite  dielectric/vacuum  interface,  this  dispersion  relation  is  found  from^^^ 


c^kii^  _ e(cl') 

e(cj)  + 1 


( B-2) 


Then  if  Op  and  the  frequency  of  the  incident  radiation  are  adjusted  so  that  for 
some  n,  u=w^(k,/'^0,  the  grating  exhibits  a resonant  response  to  the  incoming 
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field,  and  a surface  polariton  is  generated.  In  the  present  paper,  we  are 
Interested  in  inciaent  beams  that  operate  within  such  a resonance  peak. 

To  proceed,  we  introduce  a convention  that  will  be  used  throughout  the  paper. 
Its  purpose  is  to  keep  the  notation  and  analysis  as  simple  as  possible.  The  inci- 
dent electric  field  is  a vector,  but  we  work  with  only  a single  Cartesian  component. 
We  choose  the  component  parallel  to  the  surface.  Furthermore,  we  shall  need  only 
the  amplitude  of  this  component  very  slightly  above  the  surface.  The  values  of 
this  electric  field  component  elsewhere  in  space,  and  the  Cartesian  component  per- 
pendicular to  the  surface,  are  readily  constructed  from  this  information  by  elemen- 
tary considerations.  This  convention  allows  convenient  application  of  the  results 
of  Reference  (6),  since  the  discussion  there  focuses  on  the  relation  between 
electric  field  components  of  the  various  beams  parallel  to  the  surface.  Also,  this 
convention  reduces  the  present  problem  to  one  that  can  be  discussed  with  use  of 
notation  appropriate  to  wave  propagation  in  one  dimension. 

We  let  E(x,t)  be  the  amplitude  of  the  field  component  parallel  to  the  surface 
in  the  incident  wave.  Here  x is  a coordinate  parallel  to  the  surface,  with  the  x 
axis  perpendicular  to  the  grating  grooves.  We  wish  the  electric  field  at  each 
point  to  have  the  form  of  a chirped  pulse,  with  gaussian  envelope: 


E 


exp 


(jjp  t(l  + xt) 


( B-3) 


^(o)  maximum  field  generated  by  the  laser  pulse,  2c^  the  width  of  the 

gaussian  intensity  envelope  of  the  pulse,  and  x is  a measure  of  the  amount  of 

linear  chirp  applied  to  the  signal.  In  other  works^^^he  combination  Xu  is 

0 

denoted  by  r/2.  We  write  Eq.  (B-3)  in  the  form 


♦ OD 


(B-4) 
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where  after  a short  exercise  one  finds 


1 /2 

E(w)  = j O exp[  - -I  j -O*  exp[--I  (o*)^  , ( B -5a) 


wi  th 


o = o-(  1+i  2o„^Aui„) 
0 0 i 0 


-1/2 


{ B -5b) 


and  the  complex  square  root  is  chosen  so  that  Re(o)>0. 

These  expressions  give  the  electric  field  at  one  point  in  the  wave.  To  des- 
cribe the  spatial  dependence  as  well  as  the  time  dependence  of  the  incident  wave, 
the  term  which  involves  (u-w^)^  is  multiplied  by  expji  k„(aj)x  j and  that  propor- 
tional to  (w+u^)^  is  multiplied  by  expj-i  k„(u)xj  where 

kjw)  = f sin  0^  (6  -6) 


is  the  component  of  the  wave  vector  parallel  to  the  surface  for  the  Fourier  compo- 
nent with  frequency  u>. 

Finally,  to  insure  the  wave  illuminates  only  a strip  of  width  L on  the  surface, 
as  illustrated  in  Figure  (1),  we  multiply  the  final  expression  for  E(x,t)  by 
e(x+L/2)  •e(L/2-x)  , where  e(z)  Is  the  unit  step  function  that  equals  +1  for  posi- 

tive values  of  its  argument,  and  vanishes  for  negative  values. 

We  write  the  final  expression  for  thL  incident  field  in  the  form 


E(x.t)  ^E(k,,...)  e*''"''  e-'"‘ 


where  one  has 


E(k„,iu)  = i{2ir)^/2  E(^j  |o  exp| 


sin  I 

{k„-k„(->)) 

^kii""kn(uj)^ 

(8-7) 


( B -8) 


( B -8) 
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-o*  exp 


sin 

j(k„+k„  (c.) 

) 7 

kii kii  ( ui) 

The  next  task  is  to  describe  the  response  of  the  grating  to  the  incident 
field  in  Eg.  (B-7).  if  the  amplitude  u of  the  grooves  in  the  grating  is  small, 
so  the  fields  generated  by  the  incident  field  may  be  calculated  to  first  order  in  u, 
then  we  need  to  consider  the  response  of  the  grating  to  the  sinusoidal  field 
E(kii.w)  exp  [i(k„x-ut)],  and  the  form  of  the  output  signal  may  be  found  by  synthe- 
sizing the  Fourier  components  calculated  from  linear  response  theory.  This  is  the 
procedure  followed  here. 

The  interaction  of  the  grating  with  a sinusoidal  field  proportional  to 
exp[i(k„x-uit)]  may  be  constructed  from  the  prescription  given  in  Section  IV  of 
Reference  (5)  (or  from  Reference  (5)).  One  needs  the  amplitude  of  the  reflected 
part  of  the  incident  wave  and  that  transmitted  into  the  substrate.  These  are 
readily  calculated  from  the  appropriate  Fresnel  equations.  For  practical  purposes, 
we  shall  always  have  Oq“q>>1  and  Aa^«l.  Also,  the  spread  in  wave  vector  k„  impor- 
tant in  the  Fourier  integral  will  be  small  compared  to  so  long  as  L is  long 

compared  to  the  wavelength  of  the  wave  in  vacuum.  Under  these  conditions,  to 
excellent  approximation  one  can  calculate  the  reflected  and  transmitted  incident 
pulses  by  simply  multiplying  each  Fourier  component  by  the  same  Fresnel  coeffi- 
cient, namely  that  for  frequency  and  k|_=k^i(up). 

We  shall  simply  quote  the  result  for  the  form  of  the  output  pulse,  since  the 
prescription  for  computing  is  supplied  in  Section  IV  of  Reference  ( 6) . We  have 
made  some  approximations  along  the  way,  and  we  describe  them  briefly.  First  of  all, 
we  are  interested  in  frequency  regimes  where  the  dielectric  constant  e(uj)  of  the 
substrate  has  negative  real  part.  We  also  assume  Im^c(u))j<<|  Rec(w)  | , so  the  propa- 
gation length  of  the  surface  polariton  is  long.  Then  in  all  those  parts  of  the 
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response  function  which  vary  smoothly  with  frequency,  e(ui)  has  been  replaced  by 
"Up(“o)l’  magnitude  of  its  real  part  evaluated  at  This  approximation  is 

readily  avoided  if  desired,  at  the  price  of  an  increase  in  complexity  of  some  of 
the  expressions  below.  However,  for  our  purposes,  this  approximation  introduces 
little  error. 

There  is  one  factor  where  this  replacement  cannot  be  made.  In  the  calculation 
of  the  response  function,  one  encounters  the  denominator^ 


D 


(kj")  co)  = 1 + e(ui) 


a(k.,(")  w) 


( B -9) 


where  k,/'^^=k„ , d is  the  period  of  the  grating,  n the  order  of  the  dif- 
fracted beam  under  consideration,  and 


{ B -10a) 


( B -10b) 


The  quantity  o(k„^^\iu)  vanishes  when  k,/^^  = k^(iu) , with  k^{u)  the  surface  polariton 
wave  vector  (complex,  with  Re(k^(cj) ) >0)  calculated  from  the  dispersion  relation  in 
Eq.  (B-2)).  Since,  as  we  shall  see,  o(k„^'^\uj)  is  in  the  denominator  of  the  res- 
ponse function  and  we  shall  work  in  the  near  vicinity  of  this  resonance,  the  full 
complex  dielectric  constant  e{aj)  will  be  retained  here.  Also,  in  the  numerator  of 
the  expressions  for  the  diffracted  fields,  since  we  work  close  to  the  zero  in 
D(k,/^\u),  wherever  the  factor  e( uj) a^(k,/^\uj)/a(k,/^\aj)  is  encountered,  it  will 
be  replaced  by  -1 . 
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Then  if  we  presume  it  is  the  order  diffracted  beam  that  generates  the 
surface  polariton,  in  this  beam  the  tangential  component  of  the  electric  field 
E^^^(x.t)  just  above  and  below  the  surface  are  equal.  Following  Eq.(B-7)  we  write 


•foo  +0 


E<=>(x,t)  . J yf-E(^)(k,„.)  e' 


X ^-iujt 
e 


( B-ll) 


where  the  grating  response  theory  described  above  produces  the  expression^®^ 

2i  u^"^  T(u>„)  / , , 

m1  = 2.  0 r /^  2-rrn  \ 


where  D(k„,oj)  has  been  defined  above, 


T(«„) 


and 


a(o)  = l^r^^o 


)| 


( B-12) 


( B-13a) 


( B -13b) 


Finally,  if  the  function  u(x)  is  the  grating  profile,  we  write 

= nS«  ^ • (B-14) 

Thus,  the  quantity  u^^^  in  Eq,  (B-12)  is  the  amplitude  of  the  Fourier  component  of 
the  grating  profile  shape  function  u(x)  responsible  for  Bragg  diffracting  the  inci- 
dent wave,  to  generate  the  surface  polariton. 

If  the  above  results  are  collected  together,  then  after  some  algebra  which 
makes  use  of  the  identity  D(k„,  -w)  = D(k^^,u)*,  we  find 


TU.9-J 
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EW(x.t). 


X Re 


+ <*> 

0 T(to^)  ^ duj  e 


-(o) 


° I 

2 ' 0 


dk„  e " sin 


'■/- 


In  this  expression, 


(n)  _ 


. (o)  2irn 
k..  + -g- 


<*>  , ^ . 2rn 

- sin  o„  + — T- 
c 0 d 


( B-15) 


( B-16) 


The  next  task  is  to  cast  Eq.  (B-15)  into  a useful  form  that  can  be  evaluated 
for  various  cases  of  interest.  We  do  this  by  making  rather  extensive  use  of  the 
fact  that  we  are  operating  very  close  to  the  surface  polariton  resonance  in  the 
response  function  [D(k„,u)]  \ For  example,  after  some  algebra,  we  may  write  (with 
no  approximation) 


1 


D(k..,ai) 


1 


1 * f 1 “o*''-"’ 

1 + e(u; 


a(k„,u) 


( B -17a) 


«(k,.,u)  [a(k,,,to)  - e(io)  ap(k|,,u)] 

l-E^(aj) 


[k.,^  - k 2(„)] 


( B -17b) 


where 


i,  2/  \ _ u z{u) 

•'s  - ;2  [--(w)-'-“"t] 


( B-18) 


is  the  square  of  the  (complex)  surface  polariton  wave  vector,  for  a surface  polariton 
of  frequency  oj.  We  replace  k„by  kg(io)  everywhere  in  Eq.  ( B-17b)  save  for  the  factor 


i 
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*^i.  "*^s  shall  always  be  close  to  resonance.  Then  again 

everywhere  but  in  this  resonant  factor,  we  replace  c(u)  by  -1c^(uJq)|,  as  we 
have  earlier.  We  then  have  the  approximate  expression 


f 

r 


E^^^(x.t) 


le^(uQ)  |-1  ^ 


{ B -19) 


The  result  in  Eq.  ( B-19)  may  be  used  to  calculate  the  output  field  anywhere 
in  space,  provided  the  input  pulse  has  a frequency  spectrum  sufficiently  narrow 
that  the  assumptions  used  in  deriving  it  are  correct.  (A  key  assumption  is,  in 
essence,  that  u)q>>1/|oI).  The  interest  here  is  in  the  form  of  the  surface  polariton 
pulse  in  the  limit  x ->  ®,  far  from  the  region  -L/2<x<+L/2  where  the  exciting  radia- 
tion strikes  the  grating.  For  x>L/2,  the  k„  integration  in  Eq.  ( B-19)  may  be 
evaluated  by  a straight-forward  contour  integration.  This  gives  us,  again  for  x>L/2 


E^^^x,t) 


k,(.o)l 


3/2 


( B -20) 
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In  essence,  the  expression  in  Eq.  ( B -20)  is  our  principal  result.  In  it. 


k^(u)  is  that  root  of  Eq.  (lB-18)  with  lm(k^(w))>0.  We  shall  examine  the  shape  of 
the  output  pulse  through  study  of  the  structure  of  the  integrand  in  the  next  sec- 
tion of  this  paper.  We  conclude  this  section  by  casting  the  result  in  a more 
useful  form.  As  it  stands  presently,  the  integral  over  u in  Eq.  ( B-20)  would  be 
rather  tricky  to  evaluate  by  numerical  methods. 

To  make  the  notation  simpler,  define 


t(u)o)  = 


:r(^o)l'^' 

\W 


T(.^) 


{ B-21) 


and  split  the  sin  up  into  complex  exponentials  to  write  Eq.  (B-20)  in  the  form 


. {if  „(">  (^) 


X Re 


. irnL  +* 

a t(w  ) 2^  ye 

p=±l  k (u)-k„ 


where  we  have  introduced 


+ 00 

/ 


1 Z(  ,2 

do)  ■ 

[TTr  ® 


= x + y 2 


Lsine 


t^  = t + 


0 

2T~ 


X e 


i[ks(co)x^j-ujt^j] 


( B-22) 


( B -23a) 


Now  we  have  lm(k^(ui))>0  , as  remarked  earlier.  We  may  then  write 


( B-23b) 
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i / ds  e 


kj(u)-k, 


( B-24) 


SO  that,  introducing 


( B-25a) 


^ ®o 


( B-25b) 


we  have 


. unL 


Re  I i a t(w„)  S w 


° u=±l 


^ f.  r 


. 2Trns 

’ ~T~ 


( B -26) 


X / du  e 


Next  examine  the  integral  over  u.  Upon  writing  u=UQ+n  , this  integral  becomes 


I = e''“o  f dne^ 


7 -jo^  i k («  -t-n)  ^ is)  -in  t (s) 
/ dn  e ^ e ^ ° ^ e ^ 


( B-27) 


Again  with  the  presumption  that  the  incident  pulse  is  composed  of  a reasonably 
small  spectrum  of  frequencies,  \;e  expand  k^(u>Q+n)  in  a Taylor  series  about 


ks(wo+n)  = ks(uo)  + 7^  n - ]■  + - - 


where 
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(^) 

\ 3w  /u) 


( B -29a) 


3w^/u), 


( B -29b) 


In  general,  note  that  and  B are  complex  numbers,  although  they  become  real  when 
e(u)  is  real. 

After  the  expansion  displayed  in  Eq.  (B  -28)  is  used,  the  integral  over  n in 
Eq.  (B-27)  is  readily  evaluated.  We  write  the  result 


(2w)" 


? r T^(s) 

e exp = 

2o/(s) 


where 


► ^(S)  = k^(a.^)  C^(S)  -0.^  T^(S)  , 


O^^(s)  = + i B c^(s) 


T^(s)  = T^(s)-^C^(s) 


Then  after  some  algebra,  the  pulse  generated  by  the  grating  has  the  form,  with 

t(w^)  = |t(ujj)  I exp(i6)  , ( B -32] 


E^^^x,t)  = 2 E 


X - U)„  t +6 

u 0 w 


-.U.,  - 


7^  ' V*i<W 
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L>1 

2»/(s) 


( B-33) 


where  we  have 


0 • « 4.  2Tin 

— sin  0^  + — 7- 
c 0 d 


To  write  out  a^(s)  and  T^(s)  explicitly  we  have 


2 2 

a (s)  = o + i 8 X + i 6 s 

P P 


( B-34) 


( B-35a) 
( B -35b) 


It  is  possible  to  simplify  Eq.  (B-33)  considerably  through  further  algebraic 
manipulations.  To  do  this,  define 


( B-36) 


and  note  that  one 

may  write 

' 

k.(io^)  X - u„  t 
sop  Op 

3 

</) 

II 

1h 

;a. 

1 

X - lOo  ^ ^ ^ k,,(uo)  p L 

(B-37) 

We  also  have 

T^(s)  = t - 

U (s  + i p L)  , 

^ % 

(B-38a) 

and 

O^^(s)  = 0^  + 

i sjx  + s + ^ p j , 

(B-38b)  1 

where 

1 

1 

9 9 1 

- Vg  sin  e^/c| 

( B-38C)  || 
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Now  if  in  each  term  in  the  sum  over  \i  in  Eq.  ( B-33)  we  write 


C = s + y u L 


( B -39) 


then  Eq.  ( B-33)  becomes 


= 2 E(„,  u<">  (ji)  |t(..„)|  Re  i , e I = » “ I 


^ >■/ ^ 


where 


[-  ^1 ! 

( B-40) 

L 2o^(0  J i 

T(0  = t - ^ - I- 
^ ^9 


( B-41a) 


U)  = + i B[x+c]  . ( B-41b) 

In  Eq.  ( B-40),  the  sum  on  u ™y  be  carried  out  explicitly  to  yield,  after  we 


C = T L n 


Ee(>^.t)  = T E 


2 ‘^(o) 


I i(k  (o)  )x-w  t + 6| 

f Re  i 0 e • ^ ° ° I 


T^(e) 

2o^(n) 


where  we  define 


f = 2 u^"^  L 


|t(.,)l 


( B -42) 


( B -43) 


( B-44) 


. *"14  * 
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The  fom  in  Eq.  ( B-43)  should  prove  most  convenient  for  numerical  studies  of 
the  pulse  shape,  since  the  integral  that  remains  covers  only  the  finite  interval 
from  -1  to  +1.  In  Eq.  ( B-43),  we  have 


('  ■ H 


( B-45a) 


o^(n)  = 0^  + i ^ ^ Lnj 


( B-45b) 


It  is  difficult  to  obtain  a feeling  for  the  shape  of  the  output  pulse  from  the 
general  form  of  Eq.  (B-33).  '.ie  see,  however,  that  the  output  pulse  can  be  regarded 
as  a linear  superposition  of  a series  of  chirped  output  pulses.  The  question  that 
remains  is  to  decide  when  the  output  pulse  has  a shape  that  reflects  the  frequency 
composition  of  the  input  pulse,  and  how  it  is  modified  if  it  does  not.  We  address 
this  question  in  Section  C , through  analysis  of  the  general  results  here,  and 
appeal  to  a simple  limiting  case. 


-*•  7**  ••  » 
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C . The  Frequency  Spectrum  of  the  Output  Pulse;  General  Remarks 
and  a Specific  Example 

We  first  consider  the  frequency  spectrum  of  the  output  pulse,  by 
examining  the  general  structure  of  the  results  in  section  B . For  this 
discussion,  the  expression  in  Eq.  ( B -20)  will  prove  particularly  useful. 

Let  us  first  inquire  into  the  conditions  under  which  the  frequency 
spectrum  of  the  surface  polariton  pulse  launched  from  the  grating  is 
identical  to  that  of  the  input  pulse.  If  we  examine  the  structure  of 
Eq.  (B-20),  then  we  see  that  if  L is  sufficiently  small,  the  frequency 
spectrum  of  the  output  pulse  mimics  that  of  the  input  pulse.  In  particular, 
begin  by  considering  the  extreme  case  where  we  take  the  li...it  L-0.  Regard 
this  as  a formal  mathematical  exercise  at  the  moment. 

As  L — 0,  we  make  the  replacement 


sin 

( C-1) 


Then  Eq.  ( B -20)  reads 


,3/2 


[|cR(o)o)|-lj^/^  [1er(o.o)|+i] 


x Re 


t+oo 

o J"  dcije 

*00 


gi[k^(cj)x-ut] 


( C -2) 


Quite  clearly,  near  the  grating  (x  — 0),  the  frequency  spectrum  of 
the  transmitted  pulse  is  identical  to  that  of  the  input  pulse.  As  the 


'.H 
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surface  polariton  pulse  propagates  down  the  grating,  it  will  either  be  compressed 
or  it  will  expand,  depending  on  how  the  input  pulse  is  chirped,  and  the  detailed 
form  of  k^(u).  The  integral  over  w in  Eq.  ( C -2)  is  readily  evaluated  if  the 
expansion  for  k^(u)  in  Eq.  ( B -28)  is  used.  The  algebra  is  identical  to  that  in 
the  recent  study  of  pulse  compression  through  coupling  to  bulk  polari tons. ^ 

Note  that  it  is  convenient  to  break  T{uj^)  up  into  a amplitude  and  phase,  as  in 
Eq.  I B-32).  One  then  sees  that  the  phase  angle  in  T(u^)  appears  simply  as  a 
constant  phase  shift  in  the  output  wave  which  does  not  affect  the  shape  of  the 
envelope. 

The  next  step  is  to  inquire  into  the  criterion  that  must  be  satisfied  for  the 
use  of  the  limit  L -*■  0 in  the  calculation  of  the  output  pulse  shape.  This  condi- 
tion is  simply  obtained. 

We  suppose  the  incident  beam  directed  so  the  grating  is  driven  at  maximum 
efficiency.  This  requires  the  condition 


— sinG 
c 0 


to  be  satisfied,  or  that  the  quantity  Ak,,  defined  in  Eq.  ( B-36)  vanish.  We  presume 
this  to  be  the  case  in  what  follows. 

Now  the  width  of  the  pulse  in  the  frequency  domain  is  clearly  |al"\  and  we 


expect  u^|o|>>l.  Then  we  make  the  replacement  in  Eq.  ( B -20) 


(w-“o^  V~ 
9 


1)1-20 


Eq.  ( C -3).  Now  we  see  the  limit  L -*■  0 is  appropriate  so  long  as  the  inequality 


(l-  ^ sinoj  « 


( C -5) 


% 

is  satisfied.  Or  in  terms  of  the  quantity  introduced  in  Eq.  ( B-38c),  we  may 
write  the  inequality  in  the  form 


L « |olVg 


( C -6) 


Note  that  since  one  will  operate  the  grating  in  a frequency  regime  where 
is  less  than  c by  a fair  amount,  if  signal  processing  via  surface  polaritons  is  of 
interest,  we  shall  have  V =V„  in  this  instance. 

g g 

The  criterion  displayed  in  Eq.  ( C -6)  is  the  principal  limit  on  operating 
conditions  of  the  grating  coupler,  if  one  wishes  to  produce  a surface  polariton 
pulse  which  faithfully  reproduces  the  frequency  spectrum  of  the  input  pulse.  This 
limit  applies  to  a grating  operated  under  conditions  of  optimal  coupling,  in  the 
sense  that  we  have  presumed  the  criterion  in  Eq.  ( C -3)  has  been  met.  If  the 
criterion  ii  not  met  precisely,  the  conclusions  here  remain  applicable  if  the  con- 
dition 

L|ks(wo^  - I ^ ^ 


is  obeyed. 

One  can  obtain  a clear  feeling  for  the  nature  of  the  transmitted  pulse  by 
plotting  the  two  principal  factors  in  the  integrand  of  Eq.  ( B -20X  for  the  two 

*V  ^ 

cases  L<<Vg|o|  and  the  opposite  limit  L>>Vg|ol.  We  do  this  in  Figure  (2).  The 
function  sin  [(k^(w)-k,/"^)y/[k  describes  the  response  of  the  grating 

'V 

to  the  incident  beam.  As  one  can  see  from  Figure  2(a),  when  L<<|o(Vg,  one  may 


regard  the  grating  response  as  flat,  over  the  frequency  regime  appropriate  to  the 
input  pulse.  The  frequency  distribution  in  the  surface  polariton  pulse  generated 
by  the  grating  then  reflects  that  of  the  input  pulse  quite  directly.  On  the  other 

'X. 

hand,  when  L»lo|Vg,  as  one  sees  from  Figure  2(b),  the  spectral  width  of  the  pulse 
is  broad  compared  to  the  frequency  spread  of  the  incident  beam.  Then  in  this  case, 
as  is  clear  from  Figure  2(b),  only  the  central  component  of  the  input  pulse  is 
passed  by  the  grating.  Then  the  surface  polariton  pulse  generated  by  the  grating 
has  a narrower  frequency  spectrum  than  the  input  pulse;  in  coordinate  space,  the 
surface  polariton  pulse  is  longer  than  it  would  be  if  the  inequality  were  reversed 
( i f one  had  L<<| o | V^) . 

We  next  illustrate  these  points  by  considering  a simple  limiting  case  of  the 
general  formulae  above.  For  simplicity,  suppose  we  have  a gaussian  unchirped  pulse 
incident  on  the  grating.  Then  the  quantity  x in  Eq.  ( B-3)  vanishes,  while  a in 
Eq.  ( B -5b)  reduces  to  a^.  Suppose  further  we  neglect  the  quadratic  term  in  Eq. 

( B -28)  by  setting  the  coefficient  e equal  to  zero.  Then  Eq.  (6-43)  becomes 

1 i[k,(u„)x-ui„t+6] 


/.n  e’ 


n Ak„  L 


exp  - 


( C -8) 


If  we  let 


C = t - 


( C -9) 


1 I (<^n)x-u)  t+6] 

E,(x,t)  f Reje  ^ 0 “ 
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+i 

e 


( C -10) 


For  the  moment,  we  presume  the  grating  is  operated  sufficiently  close  to  opti- 
mal coupling  that  the  condition  in  Eq.  ( C -7)  is  satisfied.  In  the  notation  of 
Eq.  ( C -10),  this  becomes  Ak^^L<<l.  In  this  limit  the  factor  exp^+i  j may 

be  set  to  unity.  Later  in  the  section,  vie  investigate  the  role  of  this  factor, 
when  it  is  retained. 

The  next  step  is  to  rewrite  the  integral  over  n through  the  identity 

4*  1 OO  OD 

J"  dn  f(n)  = ^ dn  f(n)  - J"  dr\  f(n)  • ( C -11) 

-1  -1  1 

Then  the  expression  that  results  may  be  expressed  in  terms  of  the  error  func- 
tion erf(x),  defined  by 


X 

erf(x)  = J"  exp(-s^) 

^ 0 

We  find  the  result 


( C -12) 


E^'^x.t) 


(f) 


1 

E, 


sin  [k^(u,^)x-io^t+5] 


1 I I ^ ^ I 

sgn(;  t ^o^A)  erf  


( C-13) 


-sgn{ C - ^ OqA)  erf 


where  we  have  introduced 
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A = 


'X, 

°o^9 


( C -14a) 


and  we  remind  the  reader  that 


4 = t - 


( C -14b) 


The  expression  in  Eq.  ( t -13)  describes  a surface  polariton  pulse,  with  fre- 
quency and  wave  vector  The  expression  in  curly  brackets  describes  the 

envelope  of  the  signal,  which  propagates  with  the  group  velocity  V^.  The  envelope 
remains  undistorted  as  it  propagates;  this  is  a consequence  of  the  fact  that  we  set 
the  curvature  parameter  S=0.  Quite  clearly,  the  shape  of  the  envelope  is  controlled 
by  the  ratio  A=L/(a^V^) , as  our  earlier  argument  suggests.  We  explore  the  two 
limits  A«1  (where  our  earlier  argument  suggests  the  information  in  the  incident 
pulse  is  transmitted  fully  to  the  surface  polariton)  and  the  limit  A>>1  (where  the 
incident  pulse  should  be  lengthened.) 

(i)  The  limit  L<<o^V^: 

Here  the  quantity  in  curly  brackets  may  be  expanded,  with  the  first  order  term 
in  L retained.  One  finds  the  simple  result 


E^'^(x,t) 


f E 


(o) 


exp 


c2 

2a  J 


Sin  j^k^(aj^)x-u^t+6j 


( C -15) 


This  describes  a gaussian  pulse  transmitted  to  the  substrate  undistorted  by 
the  grating.  Our  approximations  in  deriving  this  result  ignore  subsequent  expan- 
sion or  contraction  of  the  pulse  by  the  curvature  in  the  surface  polariton  disper- 
sion relation.  The  gaussian  pulse  propagates  with  the  group  velocity  of  the 
surface  polariton. 


->»  f»r  . *s^.  ' V^  ••  •»  •• 
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One  by-product  of  the  present  analysis  is  an  expression  for  the  absolute 
amplitude  of  the  surface  polariton  which  propagates  on  the  free  substrate,  to  the 
right  of  the  area  illuminated  by  the  incident  beam.  Despite  the  number  of  recent 
investigations  of  the  interaction  of  radiation  with  surface  roughness  or  small 
amplitude  structures  on  the  surface,  we  have  not  seen  such  an  expression  derived. 
From  Eq.  ( C -15),  we  see  the  ratio  of  the  transmitted  surface  polariton  pulse  to 
that  of  the  incident  pulse  is  given  by  the  factor  f defined  in  Eq.  ( B-44).  It  is 
convenient  to  note  that  the  expression  for  f may  be  rewritten  to  read 


3 

/%\2  |c|^ 

1 

^|el+sin^0^  -le|^  sin0Q 

\c  / 

1 

1 

||e|-l 

P [1e|cos^0  +sin^ 

1 ' ' 

1 t ' ' J 1 ' ' 0 

ol 

( C -16) 


The  expressions  in  Eq.  ( C -15)  and  Eq.  ( C -16)  may  be  used  to  calculate  the 
absolute  efficiency  of  a grating  ruled  on  a substrate,  in  the  approximation  that  the 
amplitude  of  the  grating  is  regarded  as  small.  Note,  as  described  at  the  beginning 
of  this  discussion,  Eq.  ( C -15)  relates  the  amplitude  of  the  component  of  the  inci- 
dent field  parallel  to  the  substrate  to  the  component  of  the  electric  field 

in  the  surface  polariton  parallel  to  the  surface.  The  complete  field  of  the  surface 
polariton  can  be  constructed  by  using  this  information  in  combination  with  Maxwell's 
equations 

(ii)  The  Limit  L»a  V 
0 q 


For  general  values  of  L/o^Vg,  the  electric  field  set  up  by  a gaussian  input 
pulse  is  given  by  Eq.  ( C -13),  under  the  conditions  used  to  obtain  the  result.  To 


consider  the  limit  L>>o„V  , introduce  the  coordinate 

0 g 


X = X - V t 
9 


( C -17a) 
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and  let 


r = V / V 
9 9 


( C-17b) 


Then  except  very  near  the  ends  of  the  pulse,  when  L>>o^Vg  we  can  make  the  replace- 
ment 


1 f I ^ o"  I 

sgn(  t 1 a^A)  erf 


/2  a. 


- sgn( t 1 a^A)  erf 


' ? °o^l 
/2  a. 


( C -18) 


-V  2 rL-x  j 0^x  + ^ rlj 


The  envelope  of  the  pulse  is  now  "square".  Actually,  the  envelope  does  not 
jump  di scontinuously  zero  to  unity  at  jT=-r^  , and  from  unity  to  zero  at  x=+r^  . 
Rather,  consi derati on  of  the  full  form  in  Eq.  ( C -13)  shows  the  envelope  makes  a 


gradual  transition  from  zero  to  unity  over  a region  with  spatial  extent  Ax=a^V^<<i. 
Then  with  Eq.  ( C -18),  except  near  the  ends  of  the  envelope,  when  °QVg<<L 


we  have 


E^^^x,t) 


X sin[k^(oj^)x-u^t+()i]  . { C-19) 

We  call  attention  to  two  figures  which  serve  to  illustrate  the  results  obtained 
above.  In  Figure  3(a)  and  Figure  3(b),  we  sketch  the  shape  of  the  pulse  generated 
by  the  grating  for  the  two  cases  L<<o^Vg  and  L>>o^Vg.  In  Figure  (4),  we  plot  the 
amplitude  of  the  peak  field  in  the  transmitted  pulse  as  a function  of  L.  When 


r 
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L<OpVg,  the  amplitude  increases  linearly  with  L,  to  saturate  at  a value  independent 

of  L when  o„V„<<L. 

0 g 

We  conclude  by  exploring  the  effect  of  retaining  a finite  value  for  the  para- 
meter Ak^^L,  which  was  assumed  very  small  compared  to  unity  in  the  discussion  above. 
We  shall  be  able  to  explore  the  influence  of  this  parameter  only  for  the  limiting 
value  Ak„L»l,  but  the  results  give  an  intuitive  feeling  for  the  role  of  this  para- 
meter. Further  investigation  may  employ  the  full  form  in  Eq.  ( B -43) , with  the 
integral  evaluated  numerically.  As  before,  in  the  present  discussion,  we  assume  an 
unchirped  incident  pulse  (x=0),  and  set  6=0  also. 

We  begin  with  Eq.  ( B -43) , and  let 


(t-x/Vq) 
^ a. 


{ C -20) 


and  we  then  switch  the  integration  variable  to 


_ n . _ n L 

X - — — A - ■ ’ — 

2^  2/2’  o V 

0 g 


( c -21) 


Then  finally  we  let 


( C -22) 


to  write  the  integral  in  Eq.  ( B-43)  in  the  form 


Ej(x,t)  = /?  Ejpj  ■*  exp[ik^(U|^)x-u^t+6] 


X y dx  exp|-(x-T)^j  expjipxj 


( C-23) 
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Suppose  now  we  consider  the  limit  A<<1,  where  when  Ak  L<<1 , we  found  the 


transmitted  and  incident  pulse  had  the  same  shape.  When  A<<1,  in  Eq.  ( C -23)  the 


integration  variable  in  Eq.  ^ C -23)  is  confined  to  values  small  compared  to  unity, 


so  we  may  write 


exp  -(x-T)^j  = exp|-T^+2xT 


{ C -24) 


The  integral  is  then  readily  evaluated  to  give 


= ^(o)  ^ 


2sin hr  Ak  L 


sin[k^(iij  )x-ti)  t+6] 
s o 0 


( C -25) 


This  expression  is  identical  to  that  in  Eq.  ( C -15),  save  for  the  appearance 


unity,  if  one  works  off  the  peak  in  the  grating  response  E^  is  reduced  in  amplitude 


by  the  factor  in  curly  brackets,  but  remains  undistorted. 

Now  consider  the  opposite  limit  A>>1.  In  this  limit,  if  we  presume  we  wish  to 


consider  finite  values  oft,  then  the  upper  and  lower  limits  of  the  integration  in 


Eq.  ( C -23)  may  be  replaced  by  + “ and  - ® , respectively.  Again  the  integration 


may  be  carried  out  in  closed  form  to  give 


.(x,t)  = (2ti)^  E^pj  exp[-p^j  sin|^k^(uQ)x+2pT-u)pt+6j 


( C -26) 


Or,  to  write  out  the  explicit  form  of  the  argument  of  the  sin  function. 


Ej.(x,t)  = (2n)^E(jjj  f exp[-p^] 


r 

iin  [ ^ Ak^^ 


x-(u)  -V  Ak  )t+6 
0 g " 


( C -27) 


of  the  factor  in  curly  brackets.  Thus,  when  A«l,  and  Ak  L is  not  small  compared  to 
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This  result  may  be  compared  with  that  in  Eq.  ( C-19).  The  approximation  here 
does  not  allow  description  of  the  near-vicinity  of  the  ends  of  the  pulse.  It  is 
appropriate  only  for  the  regime  |^|<^  rL,  in  the  notation  of  Eq.  ( C -19).  We  see 
from  Eq.  ( C -19)  that  there  are  two  effects  that  arise  from  the  finiteness  of  Ak,,. 
Firstly,  the  amplitude  of  the  transmitted  pulse  is  reduced  by  the  factor  exp[-p^], 
and  its  wave  vector  and  frequency  are  shifted  away  from  ^^(uq)  and  respectively 
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D . Concluding  Remarks 

While  the  principal  conclusions  we  have  reached  are  summarized  in  the 
introductory  remarks  of  Section  A,  it  may  prove  useful  to  direct  the  reader's 
attention  to  particular  results  obtained  through  the  analysis  in  Section  B and 
Section  C . 

We  consider  conversion  of  an  incident  electromagnetic  pulse  to  a surface 
polariton  pulse  through  interaction  of  an  incident  pulse  with  a grating,  as 
illustrated  in  Figure  (1).  The  electric  field  has  the  time  profile  given  in 
Eg.  (B-3).  A convenient  and  workable  expression  for  the  form  of  the  trans- 
mitted pulse  is  provided  by  Eg.  ( B-43).  To  derive  this  expression,  we  have 
presumed  o^u^<<l , i.e.,  the  freguency  spectrum  of  the  incident  pulse  has  width 
small  compared  to  This  and  the  assumption  that  the  height  of  the  grating 
structure  is  small  are  our  principal  assumptions. 

In  Section  C , we  show  that,  if  L is  the  width  of  the  grating  exposed  to 
the  incident  beam,  then,  when  L«)a|Vg,  the  transmitted  pulse  has  a shape  and 
spectral  composition  that  accurately  reflects  that  in  the  incident  pulse.  The 
guantity  a is  defined  by  Eg.  ( B-5b),  while  is  given  in  Eg.  ( B-38b). 

A simple  expression  for  the  shape  of  the  transmitted  pulse  obtains  for  general 
values  of  |o|Vg/L,  provided  we  assume  the  grating  is  driven  close  to  the  peak  in 
its  freguency  response.  The  latter  assumption  reguires  that  Ak  L<<! , where  Ak 
is  defined  in  Eg.  { B-36). 

When  Ak„L<<l,  the  shape  of  the  transmitted  pulse  is  displayed  in  Eg.  ( C -13). 
Special  limiting  forms  are  displayed  in  Eg.  ( C -15)  (the  limit  L<<|o|Vg)  and  in 
Eg.  ( C -19)  (the  limit  L>>|o|Vg).  We  note  these  formulas  display  explicit  ex- 
pressions for  the  absolute  amplitude  of  the  transmitted  pulse.  These  results  allow 
the  efficiency  of  the  grating  to  be  calculated. 
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Finally,  when  Ak„L>>!  (the  grating  is  operated  far  from  the  peak  in  its 
response).  Eg,  ( C -25)  and  Eq.  ( C-27)  provide  expressions  for  the  trans- 
mitted fields,  in  the  limits  L<<Vglo|  and  L>>Vg|o|,  respectively. 

To  derive  the  expressions  in  Eq.  ( C -15)  and  Eq.  ( C-19),  along  with 
Eq.  ( C -25)  and  Eq.  ( C-27),  we  have  set  the  parameter  \ in  Eq.  ( B -3)  equal 
to  zero,  in  the  interest  of  simplicity.  Thus,  these  special  forms  apply  to  an 
unchirped  incident  pulse  of  gaussian  profile.  It  is  straightforward  to  generalize 
the  results  to  the  case  X4O,  if  desired. 
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FIGURE  CAPTIONS 


Figure  (IV-l)The  excitation  geometry  that  forms  the  basis  for  discussion 
in  the  paper.  The  incident  beam  is  presumed  to  illuminate  a strip  of  the 
grating  of  width  L. 

Figure  (IV-2)A  plot  of  exp[-  ^ and  sin 

j^k^(o))-k^^ , for  the  two  cases  (a)  L«Vg|o|  and  (b)  L>>V  |o|  . We 
assume  here  that  a is  real,  and  in  the  figure,  Ak^^=k^(w) -k^^ . 

Figure  (IV-3)The  shape  of  the  surface  polariton  pulse  generated  by  the 

grating  when  (a)  L<<a  V„  and  (b)  L»a^V  . 

^ ^ og  og 

Figure  (IV-4)The  value  of  the  peak  field  in  the  transmitted  pulse,  as  a 
function  of  L. 


Figure  IV-1 
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Figure  IV-2 
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V.  THEORY  OF  WAVE  PROPAGATION  IN  A MEDIUM  WITH  TIME  VARYING  DIELECTRIC  CONSTANT 


D.  L.  Mills 
Department  of  Physics 
University  of  California,  Irvine 
Irvine,  California  92717 


General  Information 

One  method  of  modifying  the  frequency  composition  of  a laser  pulse 
is  to  pass  it  through  a medium  with  a dielectric  constant  that  varies 
with  time.  The  dielectric  constant  may  be  modulated  by  subjecting  the 
medium  to  an  external  probe,  such  as  an  externally  applied  electric 
field  or  strain. 

In  these  notes  we  examine  the  theory  of  propagation  of  pulses  in 
such  a medium  and  the  interaction  of  incident  radiation  with  it.  Most 
discussions  in  the  literature  of  pulse  propagation  in  media  with  time- 
varying  index  of  refraction  are  based  on  very  simple  formulae,  and  our 
aim  here  is  to  derive  the  results  systematically.  The  main  feature  of  the 
problem  we  exploit  is  the  fact  that  in  practical  situations  the  time 
variation  of  the  dielectric  constant  is  slow  on  the  time  scale  set  by  the 
carrier  frequency  of  the  laser  pulse.  In  a typical  case  we  may  have  10.6  u 
laser  radiation  (frequency  3xlol3  Hz)  which  passes  through  a material  with 
dielectric  constant  modulated  by  an  external  field  at  megacycle  frequencies. 

For  simplicity,  we  consider  a material  that  (in  the  presence  of  the 
external  field  that  modulates  the  dielectric  constant)  is  uniaxial,  with 
the  electric  field  vector  in  the  laser  pulse  directed  along  a principal 
direction.  These  assumptions  are  made  here  for  convenience,  and  the  treat- 
ment is  readily  generalized  to  propagation  in  a medium  with  time-varying 
dielectric  constant  of  tensor  character,  with  electric  f4eld  in  the  pulse 
of  arbitrary  polarization. 
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For  the  simple  configuration  described  in  the  preceding  paragraph, 
the  electric  field  E(x,t)  and  displacement  field  D(x,t)  are  quite 
generally  linked  by  the  wave  equation 


3x^ 


3^D(x,t) 

~3? 


= 0. 


0) 


We  presume  the  magnetic  permeability  of  the  material  equals  unity,  and 
then  the  relation  in  Eq.  (1)  is  independent  of  the  manner  in  which  D(x,t) 
and  E(x,t)  are  related. 

We  assume  that  D(x,t)  and  E(x,t)  are  linked  through  the  relation 


D(x,t)  = e(t)  E(x,t) 


(2) 


where  c(t)  is  the  dielectric  constant  that  depends  on  time  by  virtue  of 
the  external  field. 

We  eliminate  E(x,t)  from  Eq.  (1)  and  find  the  equation  satisfied  by 
D(x,t).  This  leads  us  to 

isiMl  - 4il  . 0.  (3) 

c ar 

an  equation  which  admits  solutions  of  separable  form: 

D(x,t)  = e^‘‘\(t)  . (4) 


where  d|((t)  satisfies 
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9\(t) 


iity 


d,(t) 


(5) 


We  begin  by  studying  the  separable  solutions  displayed  in  Eq.  (4), 
and  then  we  examine  the  theory  of  pulse  propagation  in  the  presence  of 
the  time-varying  index. 

The  structure  of  Eq.  (5)  is  identical  to  the  time-independent  form 

of  the  Schrb’dinger  equation  which  describes  motion  of  a particle  in  a 

spatially  varying  potential.  To  see  the  correspondence,  we  replace  the 

time  variable  t in  Eq.  (5)  by  the  spatial  coordinate  x,  d|^(t)  by  f^(x), 

2 2 

the  wave  function  of  a particle  of  energy  E,  and  finally  c k /e(t)  is 
replaced  by  2m(E-V(x))/fi^. 

To  examine  the  nature  of  d|^(t)  for  the  case  where  e(t)  varies  slowly 
in  time,  we  are  led  to  pursue  a method  analogous  to  that  used  in  quantum 
mechanics  to  describe  the  motion  of  a particle  which  moves  in  a potential 
that  varies  slowly  in  space.  This  is  the  WKB  method. 

To  make  the  analogy  with  the  quantum  theory  more  direct,  split  the 
dielectric  constant  e(t)  into  a background  value  and  a time  varying 
portion  Ae(t)[e(t)=EQ+Ae(t)].  Then  let 


2 . 


\(t)  = 


Ae(t) 
+ Ae( 


m 


(6a) 

(6b) 


to  write  Eq.  (5)  in  the  form 
9^d.(t)  2 

— 2 + = 0 . (7) 

3t 


V-4 

We  want  to  solve  Eq.  (7)  under  the  assumption  that  A(t)  varies 
slowly  over  a time  interval  the  order  of  . In  essence,  the  method 
used  here  expands  the  solution  in  powers  of  (n/cjQ),  where  n is  a fre- 
quency that  characterizes  the  time  variation  of  A(t). 

To  proceed,  we  write 

d^(t)  = (8) 

where  A^,(t)  and  (f)|^(t)  are  real.  Upon  inserting  this  form  into  Eq.  (7) 
and  equating  the  real  and  imaginary  parts  separately  to  zero,  one  obtains 
two  coupled  equations  for  A,^(t)  and  (|),^(t).  These  may  be  arranged  to  read 


Suppose  that  A|^(t)  is  time-independent.  Then  Eq.  (9)  and  Eq.  (10) 
admit  solutions  of  the  form 

J e 


and 
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A|^(t)  = 1. 


(11b) 


independent  of  time. 

We  seek  the  solutions  of  Eq.  (9)  and  Eq.  (10)  when  A|^(t)  varies 
slowly  in  time.  This  solution  shall  have  the  property  that, when  A|^(t) 
becomes  time-independent,  the  solution  reduces  to  that  in  Eq.  (11). 

We  begin  by  eliminating  Aj^(t)  from  Eq.  (9).  This  may  be  done  by 
integrating  Eq.  (10): 

A/(t)<^l^(t)  = A^^(0)(^l^(0)  H (^^(0).  (12) 


In  Eq.  (12),  (ji|^  = d(j)|^/dt  and  we  normalize  the  solution  (arbitrarily)  so 
that  A|^(0)=1.  From  Eq.  (12),  we  then  have  the  (exact)  result 


(13) 


Eq.  (9)  may  then  be  written  in  terms  of  (j)|^(t)  alone: 


(<J)k) 


3 * * ^ • • • • 

2 ■^‘^k 


(14) 


The  quantity  in  square  brackets  on  the  right-hand  side  of  Eq.  (14) 
vanishes  when  Aj^(t)  is  time  independent.  This  follows  from  Eq.  (11a). 
Thus,  Eq.  (14)  provides  a means  of  solving  for  4i|^(t)  by  the  method  of 
successive  approximations,  in  the  limit  where  Aj^(t)  varies  slowly  in  time. 
We  do  this  by  writing 
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4>|(  = + • 


where</i|^^°^  satisfies 


’ (o)  \ 2 2 

>k  ) = % -Ak(t)  a 


or  4>k^°^(t)  = ck  f ^ 
0 


kf  ^ 
J n(t‘ 


where  n(t'j  is  the  index  of  refraction  of  the  material. 

Upon  inserting  the  complete  form  of  Eq.  (15)  into  Eq.  (14)  and  keep- 
ing the  terms  that  provide  the  first  correction  to  Eq.  (17),  we  have 


2 V" 'r  TiiW  [I ■ V°’ V°’]  <'8) 


or,  if  all  quantities  are  expressed  in  terms  of  n(t), 


= ‘ I ^]’ 


n 


This  correction  is  first  order  in  the  time  variation  of  n(t),  and  it 
provides  a correction  to  the  lowest  order  WKB  approximation  for  ^^(t), 
which  is  given  in  Eq.  (17).  We  shall  see  that  order 


. _ A..--  _ • 


where  n is  a typical  modulation  frequency  in  e(t),  and  the  carrier 
frequency  of  the  pulse  that  propagates  through  the  material.  We  shall 
assess  the  role  of  this  correction  in  modifying  pulse  shapes  below. 

Note  that  the  amplitude  A(t)  may  &e  found  from  Eq.  (13),  once  4)|^(t)  is 
known  to  a given  order  in  n/a)^. 

Further  corrections  to  <J)|^(t)  may  be  generated  by  carrying  out  the 
successive  approximation  scheme  to  higher  order.  In  this  fashion,  one 
may  generate  an  expansion  of  the  functions  A|^(t)  and  4>|^(t)  in  powers  of 
Note  that  we  have  yet  to  make  any  approximation  that  presumes  Aj|t) 
in  Eq.  (6b)  to  be  small  in  magnitude.  The  only  assumption  is  that  A^t) 
varies  slowly  in  time.  Furthermore,  we  have  made  no  assumption  about  the 
specific  form  of  the  time  dependence  of  n(t). 

An  Application: Pulse  Propagation  in  the  Presence  of  a Time  Varying 

Index  of  RefractioF 

We  consider  propagation  of  a pulse  which  at  time  t=0  is  tailored  to 
be  Gaussian  in  form.  At  t=0,  the  electric  field  is  given  by 


E(x,0)  = 


r 


(2tt) 


jj2  E(k)e 


(21) 


where  we  choose 


E(k)  = 


(2a) 


T72 


exp 


4? 


(22) 


SO  that  the  explicit  form  of  E(x,0)  is 


E(x,0)  = E^  exp[-ax^  + ikox] 


(23) 


i 

r 


I 


H 


I 
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Upon  noting  that  at  time  t,  we  have  E(x,t)  = D(x,t)/e(t),  the  field 
E(x,t)  in  the  pulse  may  be  written  in  the  form 


E(x,0)  = 


/ 


E(k)A^(t)e’"’'e  " 


We  examine  the  explicit  form  of  the  pulse  for  the  case  where  the 
time  varying  part  Ae(t)  of  the  dielectric  constant  is  small.  In  the 
notation  of  Figure  (6b),  we  have  Ae(t)«CQ.  For  all  times,  we  may  then 
replace  the  ratio  e(o)/e(t)  in  Eq.  (24)  by  unity.  Also,  for  small  Ae(t), 
the  amplitude  A|^(t)  may  be  approximated  by  unity.  It  is  the  time-dependent 
phase  factor  ((>|^(t)  that  plays  the  greatest  role  in  affecting  the  pulse 
shape,  in  this  limit.  We  concentrate  on  the  behavior  of  (fi|^(t)  and  its 
influence  on  the  pulse  shape.  Rather  than  work  with  e(t),  we  consider 
instead  the  time  varying  index  of  refraction  n(t),  and  write 


n(t)  = + An(t), 


An(t)«n^  always. 


For  4)|^'°^(t),  we  have 


4>|,^°^t)  = ck 


^ 

J 


ck 

I 


dt  An(t) 


If  we  presume  the  index  of  refraction  is  modulated  in  a sinusoidal 


fashion. 


■ , ~ ' ji’’'  , 
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An(t)  = An  cos(nt+(j)), 


(27) 


then 


. (o)/.N  _ ck 


If  An  is  small , in 


" 4ck 


t - ^ ^ {sin(fit+(}))-sin((j))} 

"o  “ 

the  term  in  (h)^  may  be  ignored  to  give 
[n(t)-n(o)]  = - [sin(nt+(j))-sin((())] 


(28) 


(29) 


so  that,  with  u)=ck, 

+ <}.^(^)(t)  = ^ 

We  can  see  from  this  expression  that  the  effect  of  including  the  finite 
rate  of  change  of  n(t)  (i.e.,  the  effect  of  retaining  <J>|(^^^)  is  small,  and 
leads  to  no  qualitative  change  in  the  behavior  of  the  pulse.  The  term  in 

leads  to  the  contribution  to  the  right-hand  side  of  Eq.  (30)  proportional 

to  This  has  precisely  the  same  time  variation  as  the  term  in  <ti|^^°^(t) 

proportional  to  An,  and  leads  to  no  new  effects  of  a qualitative  nature. 

2 2 2 

Under  typical  operating  conditions,  the  ratio  n^  /u  is  very  small  indeed. 

A short  pulse  which  propagates  in  a material  where  the  index  is  modulated 
slowly  may  have  a time  duration  short  compared  to  the  period  T=2Tr/n  of  the 
external  field  which  modulates  the  index  of  refraction.  Then  we  may  calcu- 
late the  form  of  the  pulse  by  expanding  the  right-hand  side  of  Eq.  (30)  for 


An  / "n  “ 

- ^ P+-S— 


'sin(i^t+d))-sin(4)) 


. (30) 
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small  times  about  t=0.  If  we  do  this  and  keep  track  of  the  influence 

2 

of  the  term  proportional  to  n , we  have 


liL 

n. 


Tl(t) 


on 


T2(t) 


(31) 


where 


T-|(t)  = t 


1 + cos(<{))  + 2^*  sin($)nt 


An 


- i ^ cos(((i)  (nt)^  + • • • 


and 


(32a) 


j sin((}>)nt  - ^ cos((j))  (nt)^ 


(32b) 


If  we  recall  that  w=ck,  then  for  E(x,t)  we  have 

v2 


E(x,t) 


2(Tra) 


w 


f 


+ 00 
dk  e 


4o 


ik 


CT, 


X-  ■ 


■'“Eir 


(33) 


We  let 
k = k^  + K 


(34a) 

(34b) 


and  presume  o/kQ«l , so  only  values  of « small  compared  to  k^  are  important. 
Physically,  this  means  in  real  space  the  packet  is  many  wavelengths  long. 

We  then  can  write 


so  that 


ik_s 


E(x,t)  = 


n n*- 

■■'■cir  ^2 
e 0 


2{TTa) 


VT 


+ 00 

dk  exp 


— 00 


X exp 


Define  the  quantity 


where  = ck^  . 


Then  upon  integrating  Eq.  (36),  and  introducing 


CT. 


1 11 
T = s + ^ A=  X - — + 4 A 


and  a defined  by  ^ = a ^ "F 

0 


we  have 


E(x.t)  = 


1/2 


Eq  exp 


[-  ofT)^] 


'*'o( 
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The  quantity  A provides  a measure  of  the  effect  of  the  finite 
frequency  n of  the  modulation  rate  of  An(t)  on  the  shape  of  the  propaga- 
ting pulse.  This  effect  shifts  the  centroid  or  "center  of  gravity"  of  the 
pulse  through  appearance  of  A in  Eq.(38).  This  effect  leaves  the  shape  of 
the  pulse  unchanged. 

The  influence  of  A on  the  shape  of  the  pulse  is  found  from  the 
appearance  of  the  factor  of  A in  the  expression  for  a , which  controls  the 
width  of  the  pulse.  We  have  for  |E(x,t)|  the  expression 


(40) 


As  the  pulse  propagates,  the  effect  of  the  presence  of  A is  to  broaden 
the  pulse.  Notice  that  A increases  linearly  with  the  time,  starting  at 
zero,  if  our  approximation  nt<<l  is  meaningful.  If  we  call  Ax(t)  width  of 
the  pulse  at  time  t,  then  Eq.(40)  leads  us  to 


Ax{t)  = 


(41) 


If  x(t)=ct/n^  is  the  distance  the  pulse  has  traveled  in  the  material, 
we  have 


A 


2 


A n cos(>i)  x(t) 


(42) 


If  the  amplitude  An  of  the  modulation  of  the  index  of  retraction 
remains  fixed,  but  the  modulation  frequency  is  increased,  the  pulse 
broadening  rapidly  becomes  more  severe;  it  is  not  possible  to  narrow  the 
pulse  by  manipulating  any  parameters  which  enter  the  present  discussion. 

Under  practical  circumstances,  the  amount  of  broadening  that  occurs 
is  miniscule,  as  some  numbers  show. 

Consider  a lOu  pulse  in  a medium  where  the  index  of  refraction  is 
modulated  at  10^^  Hertz,  a very  high  modulation  rate  indeed.  (This  is 
higher  than  one  can  use  in  practice.)  Then  = 8x10"^^.  Let  the 

pulse  be  of  0.1  nanoseconds  duration,  so  ax  = CAt  = 3 cm.  If  the  pulse 
propagates  10  cm  in  the  material,  then  (x(t)/AX(o))  'v.lO.  Finally,  for 
these  numbers  c/uqAx(o)  -x.  5x10’^  , and  f i';  incredibly  tiny. 

Our  conclusion  is  that  to  extraordinaiy  precision,  the  shape  of  the 
pulse  may  be  described  by  ignoring  and  the  tiny  pulse  broadening 

effects  incorporated  in  it.  We  need  only  retain  (}i|^^°^(t)  for  a pulse 
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with  any  envelope  shape,  to  excellent  approximation 


E(x,t) 


+ 00 


E(k)  exp 


(45) 


This  formula  reproduces  the  results  of  the  elementary  description  of  pulse 
chi rping. 


Reflection  of  Monochromatic  Incident  Radiation  From  The  Surface  of  a Material 
With  Time  Varying  Index  of  Retraction 

In  the  preceeding  two  sections  the  solutions  of  Maxwell's  equations  were 
examined  for  propagation  in  a medium  with  index  of  refraction  n(t)  that 
varies  slowly  in  time  compared  with  the  carrier  frequency  of  the  wave.  If 
we  write  n(t)  in  the  form 


n(t)  = n^  + An  cos  (nt  + <())  , 


(46) 


as  in  Eq,  (27),  and  assume  An  « n^^,  then  in  the  material  we  can  construct 
characteristic  solutions  E|^(x,t)  of  the  form 


E|^(x,t)  = exp  [i(kx  - wj^t  + X|^(t))] 


(47) 


where  uj|^  = ck/ng  and 


X|^(t)  = ^ {sin  (nt  + 4>)  - sin  {^)} 


(48) 


The  expression  in  Eq.  (47)  describes  the  space  and  time  variation  of 
the  electric  field  in  the  wave  when  n <<  oj|^,  and  when  An  <<  nQ,  two  con- 
ditions very  well  satisfied  in  the  situations  of  interest  here.  In  Section 
I and  Section  II  the  corrections  to  the  form  in  Eq.  (47)  were  obtained 
explicitly,  and  their  role  assessed,  with  the  conclusion  that  for  our 
purposes  the  effectsof  the  corrections  are  quite  insignificant. 

In  this  section  we  discuss  the  reflection  of  an  incident  monochromatic 
plane  wave  at  the  surface  of  a semi-infinite  material  with  the  time  varying 
index  in  Eq.  (47).  This  will  enable  us  to  discuss  in  Section  IV  the  nature 
of  pulse  propagation  in  such  a medium,  when  a shaped  pulse  is  incident  on 
the  surface.  The  objective  is  to  assess  whether  an  unchirped  pulse  may  be 
compressed  or  altered  in  shape  if  it  strikes  the  surface  at  some  selected 
time  tg. 

The  problem  we  face  here  is  to  devise  a manner  in  which  a monochromatic 
incident  wave  may  be  matched  across  the  boundary  to  characteristic  solutions 
of  the  form  displayed  in  Eq.  (47).  For  our  discussion  it  will  prove  extremely 
useful  to  make  a Fourier  decomposition  of  the  expression  in  Eq.  (47).  This 
may  be  done  through  use  of  the  relation 

+ 00 

exp[i  z sine]  = e^"®  J^(z)  , (49) 

n=-oo 

where  J|^(z)  = (-1)  is  the  Bessel  function  of  order  k.  Through 

application  of  this  identity  with  e = nt  + (}>,  the  characteristic  solution 
in  Eq.  (47)  may  be  decomposed  into  its  Fourier  components: 

+ 00 

E|^(x,t)  = exp[ikx  - iX|^sin((^)]  ‘ 


where  we  define 
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\ 


k 


(51) 


We  introduce  the  frequency 
compact  form 


o)|^  - nn  to  write  Eq.  (50)  in  the  more 


E.(x,t)  = exp[ikx-iX|^sin((}i)]  ^ e’"'*’  Jn^^k)  • (52) 

n=-oo 

The  magnetic  field  associated  with  the  solution  in  Eq.  (52)  will  be 
required  for  use  in  the  boundary  conditions.  The  magnetic  field  H|^(x,t) 
is  readily  computed  from  Maxwell's  Equations  through  use  of  the  v x E 
equation.  A short  calculation  gives 


H|^(x,t)  = exp[ikx  - iXksin((<,)]  ^ e"''“knt  J^(x.  ) 

n=-oo 

where 

ck 


(53) 


(54) 


We  now  turn  to  the  boundary  value  problem.  The  material  with  the  time 
varying  index  of  refraction  (Eq.  (46))  is  presumed  perpendicular  to  the  x 
axis,  semi-infinite  in  extent,  with  surface  in  the  plane  x=0.  The  material 
lies  in  the  half-space  x > 0.  The  electric  field  in  all  waves  is  parallel 
to  the  z axis. 

We  presume  the  incident  wave  has  frequency  and  wave  vector  k^  = 
where  c is  the  vacuum  velocity  of  light.  The  electric  field  in  the  incident 
wave  is 

E^^^(x.t)  = exp(ikgX  - iuQt),  (55) 

and  the  magnetic  field  is 


• I'*  - 


V-17 


exp(ikoX  - iw^t)  . (56) 

We  shall  match  this  incident  wave  to  waves  of  the  form  given  in  Eq.  (52) 
and  Eq.  (53)  across  the  boundary  surface  x=0,  through  use  of  the  tangential 
E and  tangential  H boundary  condition.  While  the  form  of  the  incident  wave 
is  at  the  control  of  the  experimenter,  so  to  speak,  in  general  we  may 
expect  the  reflected  wave  to  contain  radiation  not  only  at  the  incident 
frequency  but  also  at  the  harmonics 


(i)  = Qj  - 

on  0 


(57) 


that 

appear  in 

Eq.  (52). 

Thus 

, we  take  the 

reflected  wave 

(in  vacuum)  to 

have 

the  form 

+ 00 

E^'^^(x.t) 

n=-oo 

r(R) 

n 

exp(-ik^^x  - 

(58) 

H^'^^(x,t) 

+ 00 

= -z 

p(R) 

n 

exp(-ikonX  - 

(59) 

n=-oo 


(R)  fR) 

The  amplitudes  Eq  , E)i  , etc.  are  to  be  found  by  matching  the  incident 
and  reflected  waves  to  solutions  in  the  material  with  the  form  displayed 
in  Eqs.  (52)  and  (53). 

It  is  readily  verified  that  we  cannot  satisfy  the  boundary  conditions 
by  matching  a single  eigensolution  like  that  in  Eq.  (52)  to  the  vacuum 
waves.  Suppose  we  consider  the  wave  in  Eq.  (52)  with  the  wave  vector  k 
adjusted  so  that  as  An  ->■  0,  the  solution  reduces  to  a wave  with  frequency 
ujq  equal  to  that  of  the  incident  wave.  If  the  wave  in  Eq.  (52)  is  multiplied 
by  the  prefactor  and  with  An  0 we  match  its  Fourier  component  of 
frequency  to  those  of  the  incident  wave,  and  to  the  same  Fourier 
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coniponent  of  the  reflected  wave,  then  the  amplitude  is  determined  uniquely 
by  the  two  electromagnetic  boundary  conditions.  Then  a single  amplitude 

(D) 

' cannot  be  chosen  for  the  Fourier  component  of  the  reflected  wave  of 
frequency  n ^ 0,  which  will  satisfy  both  tangential  E and  tangential 
H boundary  conditions. 

To  proceed,  we  must  match  the  waves  in  vacuum  not  to  a single  wave  of 
the  form  given  in  Eq.  (52),  but  instead  to  a linear  combination  of  them. 

We  proceed  as  follows. 

Let  E^(x,t)  be  that  eigenwave  of  the  material  which,  as  the  time  varying 
part  An(t)  of  the  index  of  refraction  is  set  to  zero,  reduces  to  a plane 
wave  of  frequency  + mn,  where  m is  an  integer.  If  we  define 


— (jj 

c m 


o)m„  = + m«  - n«  , 

mn  0 


we  then  have,  with  , 

mn  m mn 


E;(x,t)  = e;  exp[ik^x  - ix,  sin{*)]  £ e’"^  (x.  ) , (62) 


= C exp[ik^x-iXk^sin(,^)]  £ j^^Xl  | , (63) 


and  the  electric  and  magnetic  fields  in  the  material  are  found  by  superimposing 
these  waves: 


E^(x,t)  = E E^(x,t) 

m=-oD 
+ 00 

H^(x,t)  = E H^(x,t)  . 


- * * * 
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The  coefficients  are  to  be  found  from  the  electromagnetic  boundary 
conditions.  We  turn  to  the  form  of  the  boundary  conditions. 

We  first  equate  the  amplitudes  of  the  electric  and  magnetic  fields 
which  vary  with  time  like  exp(-iwot)i  then  turn  to  the  "side  band"  fields 
with  frequency  exp(-io)Q^t) . 

Note  that  in  Eq.  (62)  and  Eq.  (63)  the  terms  with  n=m  describe  fields 
which  oscillate  like  exp(-iuQt).  Thus,  conservation  of  tangential  components 
of  E leads  to  the  requirements 

eJ°^  + E E^  exp[-im(|)-iXk^sin4,]  (Xk^,)  . (66) 

m=-oo 

and 

+ 00 

C ^m;m  • (^^) 

We  next  match  components  that  vary  in  time  like  exp[-i(uQ  - P^2)t]  to 
obtain,  since  only  the  reflected  wave  appears  in  this  channel, 

+“ 

^p*^^  = £■«  e;  exp[i(p+m)^-iXk^  sin^]  (Xk„,)  . (68) 

and 

-Ep  = e'  n„.^p  exp[i(ptm)*-Uk„  sin*]  (xp„)  . (69) 

Our  task  is  to  solve  this  infinite  set  of  equations.  In  general,  this 
seems  very  difficult.  However,  in  constructing  the  characteristic  solution 
in  Eqs.  (47)  and  (48),  we  have  presumed  that  the  frequency  n that  describes 
the  rate  of  modulation  of  the  dielectric  constant  is  very  small  compared  to 
uq.  We  can  exploit  this  assumption  again  to  simplify  Eqs.  (66)  through  (69) 
to  the  point  where  the  solution  is  evident.  Indeed,  it  is  not  worthwhile 
to  consider  any  limit  other  than  that  for  which  fj  <<  uq  without  re-examining 
the  form  of  the  characteristic  solution  in  Eq.  (47). 


I 


- - 


V-20 


When  n <<  (jjQ  in  the  quantity  defined  above  in  Eq,  (62),  we  may  with 
little  error  replace  by  “ <^o  ^ done,  = n^, 

the  index  of  refraction  of  the  undisturbed  medium.  We  also  replace  by 
X|<Q  = uq  An/(nQQ)  i Xq  to  make  X)<^  independent  of  m.  Both  of  these  approx- 
imations are  well  justified  for  the  extremely  small  values  of  JI/wq  relevant 
to  contemporary  experiments. 

These  assumptions  greatly  simplify  the  structure  of  the  boundary  con- 
dition equations.  We  introduce 

= exp[im*  - iXg  sin(*)]  (70) 

and 

exp(-ip4i)  (71) 


to  cast  the  equations  in  the  form 
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The  result  in  Eq.  (77)  is  perhaps  surprising.  Even  though  the  index 
of  refraction  of  the  substrate  varies  with  time,  so  a monochromatic  incident 
wave  excites  a non-monochromatic  disturbance  in  the  medium;  the  reflected 
wave  is  perfectly  harmonic  with  frequency  equal  that  of  the  incident  wave. 
Indeed,  we  shall  see  shortly  that  the  reflection  coefficient  has  the  value 
appropriate  to  the  undisturbed  medium  with  An  = 0,  Thus,  through  study  of  the 
reflected  radiation,  one  cannot  tell  that  the  dielectric  constant  of  the 
substrate  is  modulated  in  time.  This  conclusion  is  true  only  when  n « wq  ; 
presumably  the  amplitudes  Ep  ' are  of  order  or  smaller.  Also,  as  one 

sees  from  Section  I,  we  have  presumed  in  addition  that  An/no  is  small;  for 
example,  the  prefactor  e(0)/e(t)  in  Eq.  (24)  and  the  amplitude  Ak(t)  have 
been  set  to  unity. 

The  requirement  in  Eq.  (76)  allows  one  to  guess  the  form  of  e^.  One 
may  do  this  by  noting  the  Bessel  function  addition  theorem  (identity  9.1.75 
on  page  363  of  Abramowitz  and  Stegun) 


+ 00 


Jp  - ^o)  = E Jm  (^0^  Vp  ^ M) 

^ m=-» 


(78) 


If  we  take  the  limit  Xq  ->  X]  and  presume  that  p ^ 0,  then  we  obtain 
the  statement 


+ 00 


Jq(0)  = E Jjn  ^m+p  ° » 


m=-oo 


(79) 


upon  noting  Op(x)  vanishes  as  x*^. 

The  statement  in  Eq.  (79)  tells  us  we  may  take 


m m ' 0' 


(80) 


J 
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with  found  from  Eq.  (72)  and  Eq.  (73).  Then  Eq.  (76)  is  satisfied  iden- 
tically, and  noting  that 

E (Xo)  =1  (81) 

(n=-oo 


(use  Eq.  (79)  with  p = 0,with  Jq(0)  = 1),  Eq.  (72)  and  Eq.  (73)  become  simply 


and 


E 


0 


or 


> 

= e 


(82) 

(83) 


= 2/(1  + n^)  . 

When  the  results  above  are  assembled  the  electric  field  in  the  medium 
is  given  by  the  double  sum 


E^(x,t) 


exp[ikmx-i(m-n)((,-iu)^^t]  Jn,(Xo)  J^(Xq) 


. (84) 


Through  use  of  the  identity  in  Eq.  (49),  the  fact  that  a)^^=a)^+(m-n)n, 
and  k|^=np(a)Q+mfi)/c,  both  sums  in  Eq.  (84)  may  be  evaluated  in  closed  firm, 
to  express  E^(x,t)  in  terms  of  elementary  functions.  The  result  is 


E"(x,t)  = ^ exp[ik^jX-i(.jjt] 
o'-  -* 


exp 


i X^<sin[Qt+(|)]-sin 


[•(•■  ¥H1 


(85) 


In  what  follows,  it  will  be  convenient  to  use  the  notation 


(86) 


-WW 

ilMiiiiiiii 


4 »4«]| 


(87) 


T(x,t)  = ^ |sin|^fi^t-^j  +(f)j-sin[J^t+(fi]| 


so  that  Eq.  (85)  becomes 


2E 

E^(x.t)  = 7+^  expj^-ia)p|t+T(x,t)  - 


(88) 


This  is  the  final  result  of  the  present  section. 

As  x^o,  notice  that  x(x,t)-*-o.  Thus,  the  fields  in  the  transmitted 
wave  are  monochromatic  near  the  surface.  This  accounts  for  the  fact  that 
for  our  level  of  approximation,  the  reflected  wave  is  monochromatic. 
"Sidebands"  of  the  sort  present  in  the  characteristic  solution  of  Eq.  (47) 
require  a distance  Ax~c/u  to  build  up  to  full  amplitude.  This  suggests 
that  if  one  wishes  to  "chirp"  or  otherwise  manipulate  a pulse,  the  thick- 
ness of  the  modulator  should  not  be  smaller  than  c/Q  in  the  propagation 
direction,  if  the  incident  pulse  is  to  be  modulated  in  the  optimum  fashion. 

In  the  next  section,  we  consider  a pulse  incident  from  the  vacuum 
onto  the  substrate,  and  study  the  form  of  the  transmitted  signal  through 
use  of  the  result  in  Eq.  (88)  as  the  basic  building  block. 

Transmission  of  a Pulse  Through  a Medium  with  Time  Varying  Index 

We  wish  now  to  examine  the  propagation  of  a pulse  through  the  medium 
with  time  varying  index.  We  suppose  a pulse  of  some  prearranged  shape 
is  incident  from  the  vacuum,  to  strike  the  surface  of  the  material. 

From  the  fact  that  Maxwell's  equations  are  linear,  and  in  view  of  our 
conclusion  that  just  inside  the  surface  the  transmitted  wave  excited  by  an 
incident  monochromatic  wave  is  also  monochromatic,  we  may  solve  the  pulse 
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propagation  problem  easily.  We  simply  decompose  the  incident  pulse  at  the 
surface  into  its  frequency  components,  assign  to  each  frequency  a wave  with 
amplitude,  time  and  space  variation  given  in  Eq.  (88),  then  synthesize  the 
result  to  find  the  transmitted  pulse  profile.  We  do  this  here  and  explore 
the  behavior  of  the  transmitted  pulse  for  two  cases:  a pulse  with  delta 

function  profile  and  a Gaussian  pulse. 

Let  the  incident  pulse  at  x=0  (just  to  the  left  of  the  boundary)  have 
the  time  profile  Eo(0,t)  and  write 

+ 00 

Eo(0,t)  = f ^ . (89) 

• 00 

It  then  follows  from  Eq.  (88)  that  the  transmitted  pulse  has  the  form 

[ • 

I E"(x,t)  = J ^ Eg(u))  exp  ^ioi  jt  + T(x,t)  - ^ij  , (90) 

I.  —00  ^ 

r 

which  may  be  written  in  the  very  simple  form 

E^(x.t)  = EQ(o,t  + x(x.t)  - . (91) 

All  information  about  the  change  in  shape  of  the  pulse  is  contained  in  the 
factor  x(x,t)  in  Eq.  (91).  Consider  two  cases: 

(a)  A Delta  Function  Pulse: 

Suppose  the  incident  pulse  has  the  form  of  a delta  function  which  strikes 
the  surface  at  time  tg.  Then, 


1 

Eg(0,t) 

= Eg  6 (t-to)  , 

(92) 

SO  that 

2E„  , V \ 

i. 

E^(x,t) 

= 6 1 1 + x(x,t)  - 7 - Tq| 

' c ' 

(93) 
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The  pulse  remains  a delta  function,  and  its  position  x(t)  at  time  t 
is  found  by  solving  the  equation 

x(t)  = c(t  - to)  + c X (x(t),t)  . (94) 

Suppose  for  the  moment  that  ux/c  is  small  compared  to  unity.  Then,  in 
Eq.  (94),  x(x,t)  may  be  replaced  by 

x(x,t)  - -X  cos[fit  + 4i]  , (95) 

and  Eq.  (94)  becomes,  with  c = c/Oq 

1^1  + ^ cos[qt  + (}i]j  x(t)  = ^ (t  - to)  , 

or 

x(t)  = ^ (t  - to)  . (96) 

As  long  as  ax/c  is  small  compared  to  unity,  in  n(t)  the  time  t may  be 
replaced  by  to-  Hence,  the  delta  function  pulse  propagates  with  a velocity 
controlled  by  the  value  of  the  index  of  refraction  it  sees  when  it  strikes 
the  surface.  This  result,  which  is  an  obvious  one,  holds  only  for  values 
of  X such  that  f2x/c  <<  1.  Otherwise,  one  must  solve  Eq.  (94)  fully. 

Notice  that  the  largest  possible  value  x(x,t)  can  assume  is  2 An/no^. 

Hence,  the  greatest  amount  by  which  the  position  of  the  delta  function 
pulse  can  differ  from  the  ballistic  trajectory  c"(t  - to)  appropriate  to  the 
undisturbed  medium  is  2cAn/nofJ,  which  is  a very  small  displacement  for 
typical  values  of  An  and  q. 

(b)  A Gaussian  Pulse: 

Now  we  examine  the  behavior  of  a Gaussian  pulse  incident  on  the  medium. 
Let  the  pulse  have  carrier  frequency  uq.  width  A,  and  suppose  it  strikes  the 


m 


/•» 
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surface  at  time  to.  For  the  incident  pulse  we  have 

2 


Eo(O.t)  = Eq  exp 


(t-t  )' 

-iuo(t  - to)  - j- 

tr 


The  shape  of  the  transmitted  pulse  is  then,  from  Eq.  (90), 


E"’(x,t)  = exp  |^-iuo'l'(x,t)  - 4;^(x,t)j  , 


(95) 


(96) 


where 

ij;(x,t)  = t - to  - r + x(x,t)  . (97) 

c 

Before  we  examine  the  nature  of  the  pulse  described  by  Eq.  (96),  a word 
about  the  order  of  magnitude  of  various  quantities  is  appropriate.  The  width 
of  the  input  pulse  is  A,  and  T = 2-n/^  is  the  period  of  the  field  which 
modulates  the  index  of  refraction  of  the  material.  We  shall  always  have 
T » A,  so  x(x,t)  is  slowly  varying  over  a time  interval  the  order  of  A. 

This  allows  the  discussion  of  the  pulse  shape  to  be  presented  in  simple  terms. 

Let  us  sit  at  one  position  x»  and  examine  the  time  profile  of  the  pulse 
as  it  passes  by.  The  magnitude  of  the  electric  field  reaches  its  maximum 
at  a time  tc(x)  found  from  the  condition 


ij-(x,tc(x))  = 0 . (98) 

Since  r(x,t)  varies  slowly  over  intervals  the  order  of  A,  one  may  study 
the  time  profile  of  the  pulse  by  expanding  iii(x,t)  in  a Taylor  series  about 
tc(x): 


* 

i 


V-27 


We  use  this  expansion  and  retain  only  terms  in  the  exponent  of  Eq.  (96) 
quadratic  in  (t  - tc(x)).  Note  that 


(100) 

t&l., ' 

'tc(x) 

(101) 

we  define 

-1 

a(x)  = A 1 

(102) 

wq(x)  = 1^0  I 

(103) 

■ S (■ 

3^t\ 

^'tc(x) 

(104) 

then,  E^(x,t)  becomes 


. 2E„ 

E'(x.t)  = 


(t  - t^(x))‘ 
A^(x) 


exp  ^-iuo(x)(t  - tc(x))  -i7rr(x)(t  - tc(x))^j  . (lOJ 

Since  x(x,t)  vanishes  as  x 0,  a(x)  reduces  to  A,  r(x)  vanishes,  and  wo(x) 
reduces  to  ojq. 

The  expression  in  Eq.  (105)  describes  a Gaussian  pulse  with  width  a(x) 
different  than  that  of  the  incident  pulse.  Its  frequency  is  shifted  from 
that  of  the  incident,  and  the  presence  of  r(x)  shows  it  has  been  "chirped". 

We  find 
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3t 

3t 


^ I COS  jn  |t  - + <(ij  - COS[t2t  + i|i] 


c,-r,rnf  4.  x-\ 

" sinLiJt  ■*■  i^J  _ > 

0 c 

where  the  last  statement  is  valid  when  M/c  <<  1.  Also, 


3^t 


< sin[nt  + ()>]  - sin 

n /t  - 

( 

\ c/ 

Ann 

n. 


“ ^ COS^i^t  + _ , 

0 c 


(106a) 

(106b) 


(107a) 

(107b) 


where  again  the  second  statement  is  »alid  for  nx/c. 

These  results  show  that  the  width  of  the  pulse  is  given  by,  for  nx/e  < 1, 


‘W  = tj  • (’08, 

so  the  pulse  may  be  either  broadened  or  narrowed,  depending  on  the  sign  of  the 

time  derivative  of  the  index  of  refraction  when  the  pulse  strikes  the  surface. 

When  (3n/3t)^  is  negative,  the  pulse  narrows.  The  trailing  edge  of  the 
i-o 

pulse  strikes  the  medium  at  a time  for  which  the  index  is  smaller  than  when 
the  leading  edge  strikes  the  surface.  Thus,  the  trailing  edge  "catches  up" 
with  the  leading  edge. 

The  result  in  Eq.  (108)  applies  only  when  nx/c  is  small  compared  to 
unity.  More  generally,  we  have 


a(x)  = A[1  + f(x)] 


■1 


(109) 


From  Eq.  (106)  it  is  evident  that  f(x)  does  not  increase  indefinitely,  but 
rather  its  maximum  possible  value  is  2An/nQ.  Since  An/no  is  small  compared 
to  unity,  the  fractional  amount  by  which  a pulse  can  be  compressed  by  this 
mechanism  is  quite  small. 


I 
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For  the  chirp  parameter  r(x),  again  when  ,^x/c  is  small  compared  to 
unity,  one  obtains 

r(x)  = 2^  cos[nt  + <^]  ^ . (110) 

0 c 

with  r(x)  assuming  a maximum  possible  value  A n/irn^.  The  rate  of  chirp 
r(x)  rises  linearly  with  x near  the  surface,  as  one  sees  from  Eq.  (110),  to 
assume  its  maximum  value  at  a distance  the  order  of  c/o.  from  the  surface. 
Thus,  the  ability  to  "chirp"  a pulse  will  be  inhibited  by  use  of  a sample 
that  is  too  thin. 


» • I*-** 
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ABSTRACT 

We  discuss  the  theory  of  second  harmonic  generation  by  surface 
polaritons  which  propagate  on  the  surface  of  a nearly  free  electron 
metal,  with  emphasis  on  the  infrared  frequency  range. 

resume" 

Nous  discutons  la  theorie  de  production  de  la  deuxieme  harmonique 
par  polaritons  de  surface  sur  les  metaux  dans  la  theorie  des  electrons 
presque  libres,  en  mettant  1 accent  sur  les  frequences  infra-rouge. 
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I.  Introductory  Remarks 

While  there  has  been  considerable  interest  in  the  study  of  surface 
polariton  propagation  on  a variety  of  materials^^ ^ , rather  little  attention 
has  been  devoted  to  nonlinear  interactions  between  these  waves. 

One  does  not  expect  nonlinear  effects  to  be  large,  unless  the  process 

in  question  proceeds  under  conditions  of  near  phase  matching.  One  example 

of  such  a process  has  been  proposed  by  Maddox  and  Mills^^^,  in  a study  of 

nonlinear  interactions  between  surface  polaritons  on  the  surface  of  doped 

(n  type)  zinc  blend  semiconductors.  For  a range  of  carrier  concentrations 

with  plasma  frequency  below  that  of  the  LO  phonon,  one  may  mix  two  waves 

on  the  lower  branch  of  the  two  branch  dispersion  relation  to  produce  an 

output  wave  very  near  a point  on  the  upper  branch.  This  process  has  sub- 

(3) 

sequently  been  reinvestigated  by  Bonsall  and  Maradudin'  \ who  find  an 
appreciable  output. 

In  this  note,  we  discuss  a second  nonlinear  process  which  proceeds 
under  conditions  of  near  phase  matching.  This  is  second  harmonic  generation 
by  surface  polaritons  that  propagate  on  metal  surfaces.  For  frequencies 
well  below  the  plasma  frequency  Up  of  the  conduction  electrons,  the  disper- 
sion relation  is  nearly  linear.  Thus,  the  driven  output  wave  of  frequency 
2oj  and  wave  vector  2k^^,  where  w and  are  the  frequency  and  wave  vector 
(parallel  to  the  surface)  of  the  initial  wave,  has  a wave  vector  very  close 
to  that  of  the  freely  propagating  surface  polariton,  of  frequency  2w. 

However,  it  is  unclear  at  the  outset  whether  or  not  one  can  expect 
appreciable  second  harmonic  generation  for  this  case.  At  these  low  frequencies, 
the  fields  in  the  wave  extend  quite  far  into  the  vacuum,  while  in  the  medium 
where  the  nonlinear  mixing  occurs,  the  fields  are  nonzero  only  within  the 


small  skin  depth;  nearly  all  the  wave  energy  is  stored  in  the  vacuum,  where 
no  harmonic  generation  occurs.  On  the  other  hand,  for  driving  fields  of 


given  magnitude,  the  nonlinear  currents  within  the  skin  depth  and  at  the 

(51 

metallic  surface  increase  dramatically  as  the  frequency  is  lowered.'  ' 

The  present  study  has  been  motivated  by  the  above  remarks,  and  by  an  experi- 
mental search  for  the  second  harmonic  of  a surface  polariton  propagating 
on  copper  metal,  at  the  CO2  laser  frequency. 

1 1 . Theoretical  Derivation  of  the  Second  Harmonic  Intensity 

(5 1 

In  a recent  paper,  Rudnick  and  Steam'  ' have  presented  a phenomenolo- 
gical description  of  the  second  harmonic  currents  generated  at  the  surface 
of  a simple  metal,  and  within  its  skin  depth,  when  the  metal  is  irradiated 
with  an  electromagnetic  wave  of  frequency  u.  Their  prescription  forms  the 
basis  of  the  present  study.  We  note  that  at  the  infrared  frequencies  of 
interest  here,  interband  processes  should  make  a negligible  contribution 
to  the  second  harmonic  intensity. 

We  presume  a plane,  semi-infinite  metal  in  the  half  space  z < 0.  The 
incident  surface  polariton  (frequency  u,  wave  vector  k,,)  propagates  parallel 


ir 


2 2 

e - w /id  , which  is  real 

00  P 


to  the  X axis,  with  fields  E;  ^ (z)  exp(ik /.x-iut).  The  frequency  dependent 
dielectric  constant  of  the  model  metal  is  e(u) 
and  presumed  negative  at  both  id,  and  also  2id.  Indeed,  in  the  infrared  frequency 
regime,  e((d)s -iDp/u  and  |e(id)1  » 1,  a condition'  ' we  use  below  to  simplify 
the  final  result  substantially. 

Rudnick  and  Steam  describe  three  contributions  to  the  second  harmonic 
current  density  These  are 

(i)  a surface  current  normal  to  the  surface  which  we  write 

*5(2+)  (0-))^ 


ISirmu)'' 


(1) 
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where  a is  a dimensionless  constant  of  order  unity,  m the  conduction  electron 
mass,  5(z+)  = 6(z-e)  where  e is  a positive  infinitesimal,  and  (0-)  the 
z component  of  the  incident  field  just  below  the  surface  of  the  metal. 

(ii)  a surface  current  parallel  to  the  surface  given  by 


i(2oo) 

'^SX 


(x) 


Airmu'' 


e(“)(o-)  E^“^0-)  6(z+)  (2) 


where  b is  a dimensionless  constant  of  order  unity,  and 

(iii)  a volume  current  that  exists  throughout  the  skin  depth 


leirmu)-^  ' 


-i2ut 
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(3) 


To  construct  the  fields  associated  with  the  second  harmonic,  we  insert 
the  total  current  density  given  by  Eqs.  (1)  - (3)  into  Maxwell's  equation, 
and  use  the  Green's  function  method^®^  to  find  the  electric  field  in  the 
wave.  Since  the  Green's  function  method  has  been  employed  in  several  recent 
calculations,  in  this  note  we  omit  a description  of  the  details  and  present 
only  the  final  result. 

In  the  vacuum  above  the  crystal,  there  is  a second  harmonic  field 
gi ven  by 

E^^“^(x)  = |e(2o))  exp[i2k^^x-i2wt-ajj(2o))z]  (4) 


with  e[^“^x)  . 2ik^^  E^^‘^VaQ(2aj)  , where  a^j(2a))  = [4k^^  - 4w^/c^]  ~ 

2u)/  |e(2(jj)l^c,  with  the  approximate  form  of  valid  when  |e(2uj)|  » 1. 

The  amplitude  e(2w)  is  given  by,  in  the  limit  appropriate  to  the  present 
discussion, 
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g(2u) 
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e .p  6 
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where  6(w)=  c/Wp  is  the  skin  depth  for  the  incident  surface  polariton, 

a(w)  = a(2u))  + e(2o))  a^(2cj),  with  a^(2aj)  defined  above  and  a(2w)  = 

2 2 2^ 

[4k^/-4£(2u))u)^/c'^]^,  Eq  is  the  magnitude  of  the  electric  field  component  of 
the  incident  surface  polariton  parallel  to  the  surface,  evaluated  at  the 
surface,  and  finally  y(w)  is  a dimensionless  number  given  by 

‘ (m^f  - i (mrf  ■ 

In  the  limit  of  interest  here,  (|e(w)|  and  | c(2oj)  | >>1 ) , one  may  show  that  a((j) 

is  well  approximated  by  a(u)  s |£(2w)|'f  /L  (u),  where  L (u)  = [2k  (o))-k  (2oj)]'^ 

^ ^ till 

is  a coherence  length  that  controls  the  strength  of  the  non-linear  interaction. 

In  the  formula  for  L (u),  k..((ij)is  the  wave  vector  of  a surface  polariton  of 

^ II 

frequency  to  (this  equals  that  of  the  incident  surface  polariton),  while  k^^(2co) 

is  that  of  a surface  polariton  of  frequency  2to. 

From  the  form  of  y(w),  with  the  constants  a and  b of  order  unity,  one  sees 

the  volume  currents  contribute  to  the  non-linear  mixing  (third  term  of  Eq.  (7)) 

on  equal  footing  with  the  surface  currents.  This  is  in  contrast  to  second 

harmonic  generation  from  plane  waves  incident  on  the  boundary,  where  the  con- 
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tribution  from  the  volume  currents  is  small.'  ' Also,  note  for  a perfectly 
phase-matched  interaction,  2ky^((jj)=k^^(2o3)  and  L^(uj)  becomes  infinite.  In  our 
derivation,  we  have  implicitly  assumed  the  distance  of  travel  L of  the  incident 
wave  on  the  surface  is  long  compared  to  L^(o)).  In  the  limit  L«Lj.(aj),  in  the 
formulae,  L^(u)  should  be  replaced  by 
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In  addition  to  the  fields  in  the  vacuum  above  the  metal,  there  are 


non-zero  second  harmonic  fields  within  the  metal.  Here  one  has  a driven 


2 2 2 i 

wave  (fields  proportional  to  exp[i2k^^x+2a(uj)z]  where  a(a))  = [k^^  -w  e(u))/c  y) 


and  a wave  with  amplitude  proportional  to  exp[i2k^^x+a(2co)z] . These  fields 


are  non-zero  only  within  the  skin  depth,  and  are  similar  in  magnitude  to  the 


fields  in  the  vacuum,  which  extend  into  the  vacuum  a distance  very  much  larger 


than  the  skin  depth.  The  fields  in  the  skin  depth  thus  carry  only  a very 


small  fraction  of  the  energy  in  the  second  harmonic  fields,  and  we  ignore 


their  contribution  to  the  energy  flow  calculated  below. 


We  presume  that  the  incident  surface  polariton  propagates  along  the  sur- 


face in  a strip  of  width  d,  and  the  second  harmonic  fields  are  confined  to  a 


strip  of  the  same  width.  If  I(2u))  energy  per  unit  time  which  flows 


along  the  surface  in  the  second  harmonic,  and  I,  \ that  in  the  incident  wave. 


we  find 


64tt  e^ 
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/ a.(2o))^  rMuOr 


In  our  limit  le((jj)|  » 1 and  |e(2o))|  » 1,  this  result  becomes 


T _ p ■ c'~'  I ^ m2  T 2 


which  is  the  final  result  of  the  present  section. 


We  remind  the  reader  that  to  obtain  Eq.  (8),  we  have  begun  with  a 


general  expression  for  the  output  fields,  and  utilized  a series  of  approxi- 


mations valid  when  both  |e(o))|  and  |e(2o))l  are  large  compared  to  unity.  The 


■ • 
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input  wave  and  output  radiation  are  associated  with  the  nearly  linear  portion 
of  the  surface  polariton  dispersion  curve.  Our  presumption  is  that  only  this 
frequency  regime  is  of  potential  interest,  since  only  there  is  the  inter- 
action nearly  phase  matched, 

III.  Concluding  Remarks  and  Numerical'  Estimates 

2 2 

We  first  simplify  Eq.  (8)  further  by  presuming  that  e(o))  = e^-Wp/w  , 

2 2 

and  to  « (Op  so  in  fact  c(tj)~-tOp/to  . For  the  skin  depth,  one  has  6((o)^c/tOp 

2 3 

independent  of  frequency,  while  the  coherence  length  L ((o)scco  /3to  . Finally, 

c p 

2 2 
|y(oo)|  approaches  the  frequency  independent  limit  |y1  . where  y = (4b-a-2)/2. 

Eq.  (8)  then  reduces  to 


64tt  e^to^ 

(2(o)  ■ 9d  2 4 3 
m c (Op 


(9) 


The  result  in  Eq.  (9)  shows  that,  as  to  decreases,  for  fixed  input  power, 

the  energy  per  unit  time  carried  by  the  second  harmonic  decreases  as  (o  . The 

coherence  length  L (to)  does  not  increase  fast  enough  as  to  decreases  to  over- 

i o 

whelm  the  rapidly  increasing  factors  |c(2to)|  and  (e(to)|  ''n  ^9-  (8).  These 

two  factors  have  their  origin  in  the  fact  that,  as  to  decreases,  a progressively 
smaller  fraction  of  the  energy  in  both  the  driving  wave  and  the  second  harmonic 
are  stored  in  the  medium  where  the  non-linear  mixing  occurs. 

For  Cu  metal,  one  has  TitjOp^lOeV,  and  for  a surface  polariton  at  the 

_2 

COp  laser  frequency,  rito=0.1eV.  Thus,  to/to  slO  and  the  coherence  length 
L^(to)  ^0.5cm.  This  means  the  standard  two-prisms  ATR  geometry'  \ with 
prisms  about  0.5cm  apart,  should  prove  close  to  an  optimum  geometry  to  ob- 
serve the  second  harmonic;  the  attenuation  length  of  the  waves  is  about  0.5cm 
in  practice,  so  spacing  the  prisms  this  distance  apart  also  places  them  about 
L^((ij)  apart. 
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If  I(2(jj)  expressed  in  watts,  and  the  beam  width 

d is  in  cm,  Eq.  (9)  gives  the  quantitative  estimate 


T „ 3x10"^^ 
M2o))  ■ “d~ 


(10) 


To  see  whether  the  second  harmonic  is  detectable  in  the  infrared, 

consider  the  numbers  appropriate  to  McMullen's  attempt  to  observe  the 

second  harmonic.  The  CO2  laser  employed  here  had  an  irradiance  of  50 
2 

megawatts/cm  , in  a spot  diameter  of  100  microns.  Thus,  we  take 

_? 

d = 100  microns  = 10  cm,  while  the  total  energy  in  the  incident  beam 

is  the  irradiance  multiplied  by  the  spot  area;  after  the  beam  reflects 

off  the  coupling  prism,  about  10%  of  this  energy  is  carried  by  the 

surface  polariton.  Thus,  we  estimate  I^^j  = 380  watts,  and  Eq.  (10) 

gives  l(2o))  ~ (5^10”  ) ItI  watts.  If  10%  of  the  energy  in  the  second 

harmonic  exits  from  the  output  prism,  one  should  find  (0.5x10  ) |y| 

watts  at  the  detector.  The  limit  of  detectability  in  this  experiment 
-12 

was  about  10  watts. 

While  the  numerical  value  of  y is  unknown,  but  presumably  the  order 
of  unity,  the  present  calculation  suggests  that  the  results  of  the 
experiment  outlined  above  were  inconclusive;  with  a bit  of  luck,  the 
second  harmonic  might  have  been  detected.  The  theory  shows  the  second 
harmonic  output  to  be  modest,  but  of  detectable  magnitude. 
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